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NOMENCLATURE 


Surface area of heater 
Aspect ratio, lAi 

Spacing between tvro parallel plates 
Acceleration due to aravity 

3 

Grashoff number, — 

Protrusion of heater 
Heat transfer coefficients 
Themal conductivity fluid 

Characteristic length of heater (vertical length 
parallel to flow) 

Length of test cell parallel to the test beam 

Nusselt number, hi/k 

Prandtl number, pCp/k 

Total heat rate 

Heat flux 

Resistance of heater 


Rayleigh number, . 

Reynolds number, ^ 
facing between the two heaters 
Tenperature of heated chip surface 
Ambient temperature 
Mean tenperature 

Centre line velocity 
Velocity in axial direction 
Voltage 

Velocity in transverse direction 



X 


Axial coordinate 


y = Transverse coordinate 

^ = Kinematic viscosity 

a = Thermal diffusivity 

C = Specific heat 

P 

P = Density 

M- = Viscosity 

P = Thermal es^ansion coefficient / K 

c 

£(::) = Fringe width at x location 

^ = Emissivity 


DUbscrir 


1 = heater It upper heater 

I = heater 2, central heater 

} = heater 3# loxN^er heater 

::onf iqurat! ons (Multichip heater) 


configuration refers to 

C-H-C 
C-H-H 
H-H-K 
H-H-C 
H-C-H 

.11 the figures have dimensions in mm® 


heater 1 

heater 2 

heater 3 

cold 

hot 

cold 

cold 

hot 

hot 

hot 

hot 

hot 

hot 

hot 

cold 

hot 

cold 

hot 
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abstract 

The present work aims at determining the flow patterns and 
Tiperature field in arrays of heat generating circuit boards, Che 
rcuit boards are modelled experimentally and the relationship 
etween maximum terrperature and regime of flow (such as free, 
xed or forced convection) are sought to be determined. The 
feet of wall protrusion on flow and heat transfer are studied, 

1 the present work interferometry is used to get qualitative 
iformation about the flow field and to a cross check the results 
obtained by conventional techniques. In particular visualization 
: isotherms in the space between vertical flat surfaces with a distri- 
ition of discrete wall heaters using a Mach-Zhmder interferometer 
s one of the objectives of the present study. 
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Chapter 1 
INTRODUCTION 


The cooling of electronic equipment conbinuss to be 
an actxve field of investigations in heat transfer. Effe- 
ctive cooling of electronic equipment has become increasingly 
important as power dissipation and component density continue 
to increase substantially with the fast growth of electronics 
technology. This importance arises from reliability consi- 
derations/ the desire to decrease the transit time of ele- 
ctronic signals and the advent of new computer configurations 
and measurement devices. Reliability revolves about control 
of the temperature levels of the critical components/ while 
the transit time is diminished by shorter transmission dis- 
tances/ giving rise to closely positioned components and high 
densities of dissipated electrical power. 

As solid state circuits are miniaturized large amount 
of thermal energy is generated due to flow of current-through 
gates and connections of smaller areas within a silicon chip. 
As computing and gate switching speed increases in a chip/ 
heat generation also increases. Circuit miniaturization has 
advanced to such an extent that heat dissipation capability of 
system design has become one of the primary limiting factors. 
The main consideration in thermal design is to minimize the 
maximum temperature at the chip for a given set of design 
criteria. Recent advances in computing speed and improved 



packaging density in modern digital computers/ can be attri- 
buted to convection heat transfer research resulting in 
novel thermal design for handling high heat dissipation coeffi- 
cients. The life of an integrated circuit device and that 
of other electronic components can be prolonged if the oper- 
ating temperature is maintained at a value lower than the 
maximum specified by the manufacturer. As an estimate long- 
term reliabil Ity can be increased up to 505^ for each 18^f 
(10*^0) that a device is maintained below its maximum oper- 
ating temperature. 

To give an example of how large the heat flux can be 

consider the following numbers quoted by Oktay et al (l986). 

A power dissipation of lOW on a 5x5 mm^chip yields a heat flux 

of 5, EOS w/m^ which is only 20 times less than that on the 

o 

surface of sun. The sun's surface temperature is 6000 C while 

o 

the design temperature of today's chips is 100 C. 

Theoretically, if a chip working at this heat flux 
could somehow be tethered in still air at room temperature, 
its temperature would exceed 6000*^C, assuming only natural 
convection and no radiation, with forced convection, the 
chip can be cooled to lOOO^C with common gases and to about 
50°C with common liquids. If on the other hand, the heat is 
conducted away by means of a copper stud pressed against 
the chip, the chip's temperature can be lowered to about 30°C, 
assuming that the copper stud is held at 25°C and that the 
contact between stud and chip is perfect. If a small air gap 
of say, 1 ram is introduced between the same chip and stud, the 


Immersion cooling 


chip will reach a temperature of 500°c. 
with phase change can lower the chips temperature to approxi- 
mately 30*^C if water is the coolant, or to 45°G if fluoro- 
carbons are used. 

Thermal engineers face the challenge of removing high 
heat flxixes by means of any of the heat transfer modes- 
natural or forced convection, conduction, or phase change - 
while maintaining relatively low component temperatures. The 
working temperatures of electronic components are dictated 
not only by their inherent performance requirements but also 
by stringent reliability demands. 

Conventional Methods of Cooling 

Four main considerations while designing any cooling 
system for electronic components are: 

1. To maintain device temperature at low enough level 
to allow the system to operate. 

2. To meet reliability goals for a wide range of environ- 
ment conditions, 

3. To hold the temperature variations within a fixed 
limit . 

4. To meet requisite cooling requirements while ful- 
filling stringent acoustic noise restrictions. 



On this basis we can have following methais for the cooling 
of an electronic equipment. 

1, Natural Convection Cooling t This method is both silent 
and very reliable. It is cheap and does not introduce any 
kind of electromagnetic interference with the operating 
signals. Natural convection is the least effective way 

of removing heat. Moreover design considerations often 
mitigate against effective natural convection. Equipment for 
home and office is often laid out in compact horizontal 
packages. Open construction with vertically oriented heat 
generating components is best suited for natural convection. 
Prediction of natural convection air flow and cooling rates 
with less than ideal geometries and restricted flow paths 
poses a major challenge. In practice, circuit boards are 
separated by supporting shelves. Energy released from 
electronic chips is dissipated directly to the fluid by 
convection and to the wall by conduction. Heat conducted 
to the wall is eventually dissipated by convection to the 
adjacent channel, 

2. Forced Cooling : At places with space restriction and 

high heat flux levels, forced convection is used to achieve 
acceptable cooling performance. A fan or a blower is 
used to blow the air or cooling fluid past the heated chips 



mounted on the circuit boards. It results in a high value 
of the heat transfer coefficient and better cooling. Because 
of different geometries and flow paths in different systems 
the values of various flow parameters and pressure drop can 
not be calculated readily from first principles. Forced 
cooling is used in many personal computers and nearly all 
small computers. In a system box of a typical personal 
computer, space is at a premium and the airflow path is 
circuitous. Ccmponents that are sensitive to heat are 
located in several areas of the box, all of which must be 
reached by adequate air flow. 

Besides giving better cooling and sufficient reliabi- 
lity this method of cooling has got some limitations too. 

The fan installed for cooling occupies sufficient space. It 
increases the cost of the system as the initial cost and 
running cost of the fan has to be taken into consideration. 
The fan also produces its own electranagnetic field which 
may in some cases, lead to the corruption of signal being 
received or transferred by the electronic equipment. A 
large capacity fan also produces objectionable noise. 

Even with these limitation forced cooling is being 
adopted rapidly as the size of the IC chip is reducing day 
by day raising the heat flux level in electronic equipment. 

III. Cooling with Fin Arrays : Pin arrays are widely used 
as heat sinks for cooling electronic equipment in natural 



forced convection cooling. The most common configurations 
involve vertical surfaces having fins attached to them. The 



Fig, 1.1 : Cooling with fin arrays 

heat is transferred to these fins first and from there is 
transferred to the air passing over it. Because of the 
increase in surface area better heat transfer is achieved. 
The heat transfer enhancement is also due to the disruption 
of thermal boundary layer development. Heat transfer from 
electronic devices may be enhanced by inserting ribs in 
the flow increasing convection heat transfer from downstream 
modules . 

IV. Liquid Cooling : 

(a) Indirect Liquid Cooling ; This may also be termed 


as conduction coupled liquid cooling. The cooling of high 



performance electronics has been complicated further by 
increasingly close placement of components having high 
heat fluxes as in latest multichip modules. For the 
effective cooling these thermal resistances must be 



Fig. 1.2 ; Liquid cooling module 

reduced from the internal heat source to external final 
heat sink. One of the methods of indirect liquid cooling 
is that in which the cold plate and the chip are fabri- 
cated as an integral unit. The structure resembles a 
miniature# longtidunally finned heat sink. Closure of 
the channel is achieved with a cover plate, and liquid is 
continuously pumped through the channel. 

(b) Jet Impingement Cooling ? The specific consideration 
of impinging liquid jets at hot spots to electronics cooling 
is relatively recent. An extensive study of single phase 
jet impingement cooling for microelectronic applications has 
been done by jij and Dagan (1987) . It reduces the thermal 



resistance appreciably resulting in high values of heat 
transfer coef f icientso 

(V) Heat Transfer with Phase Change ; Liquid immersion 
cooling in a pool or forced convection boiling is another 
method of cooling of very large systems having relatively 
high heat dissipation. In this method a chemically inert 
coolant of high dielectric strength and low boiling point 
is required which limits its applications. Boiling in 
dielectric fluids is also characterised by small values of 
critical heat flux and additional problem could include 
cyclic thermal stresses and electrical noise related to 
temperature fluctuations, device contamination due to 
impurities and increased comple^iity associated with the 
need of sealed chambers. In a two phase system a thermal 
stock can arise due to a sudden drop in surface temper- 
ature that occurs upon the incipient boiling. This pheno- 
menon, called hysteresis, was encountered by Samant and 
Simon (l968) in experiments with flow boiling. 

5'orced convection boiling may occur when cooling is 
effected by an impinging jet, a channel flow, or a falling 
film. Impinging jets may be directed normal to the surface 
or at any angle off the normal. Forced convection boiling 
may also occur in free falling liquid films used to cool 
vertically mounted electronic chips. 



Two Phase Thermosyphons ; Another kind of two phase liquid 
immersion cooling is complete encapsulation of the electro- 
nic assembly. This method was first used for high power 
components and air borne electronic packages and was subse- 
quently advocated for densely populated PCBs and chip arrays. 
Saturated or subcooled boiling occurs at the surface of the 
electronic component and condensation takes place within 
the liquid or at an encapsulated cold plate. 

Chips are mounted to a substrate forming one wall 
of a module filled with a fluorocarbon. Heat is transferred 
from the chips to the coolant and from the coolant to the 
internal fins of a water cooled cold plate or an air cooled 
heat sink. Use of submerged condenser facilitates mainte- 
nance of subcooled conditions in the coolant. 



Silicon 

Chip(s) 


Module 



Pig, 1,3 : Encapsulated two phase thermosyphons 


The future electronic systems will have higher heat 


fluxes at the chip and packages level and cooling will be a 


challenge. Although no new theories are required but 
basic work on some of the heat transfer mechanisms is 
needed. Boiling phenomenon on highly polished surfaces 
such as silicon chips ought to be better understood. More 
work is needed on condensation and cryogenic heat transfer 
as they apply to electronics cooling. The importance of 
thermal interfaces and contact resistances calls for an 
understanding of non elastic contact regimes as a function 
of pressure. Better methods of analysis and simulation 
are also required. As power dissipation grows, use of 
extended surfaces such as heat sinks will become more wide 
spread. A significant progress has been made in establish- 
ing a convection heat transfer knowledge base pertinent to 
problems of electronic cooling. It covers a wide range of 
geometries, with such results it is possible for a packag- 
ing engineer to assess, at least in an approximate sense, 
the relative merits of various options for a particular cooling 
problem. However the knowledge in this field is far from 
complete. Research needs are extensive in the field of air 
cooling systems through the development of advanced heat 
sink and/or flow distribution schemes. Significant impro- 
vements may also be made in the current, state of art 
indirect cooling systems. In direct cooling systems, 
literature is sparse and research needs to encompass a broad 
spectrum. Single and two phase thermosyphons offer excellent 

opportunities for achieving moderate to large heat fluxes 
under passive operating conditions. Research is needed to under- 
stand related flow and heat transfer effects better and to 
identify optimal design and operating conditions* Moderate to 



large heat fluxes may be achieved under Icnv flow mixed convection 
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conditions, v/hile vary large fluxes (200 M/cm ) require single pha; 
forced convection or forced convection boiling. In each 
case an improved understanding of related flow phenomena 
is needed for many possible geometries. Extended surface 
arrangements that offer significant heat transfer enhance- 
ment should be identified. Emphasis should be given on the 
schemes that minimize adverse effect associated with incre- 
ased pressure losses and diminished reliability. There is 
also a need for an improvement in understanding of effects 
related to the placement of discrete heat source on a subs- 
trate. In many cases the flows are three dimensional and 
conjugate effects are significant. There clearly is a 
need for experimental as well as numerical work to be done 
in this field applied to various problems. The general 
analysis programs that are available are complex and cumber- 
some. Relatively easy, interactive programs are needed not 
only for conduction but also for convection and boiling. 

Accurate airflow analysis codes for both free and forced 
con-vection must be able to handle the complex conditions 
encountered in electronics. Greatly expanded data bases 
and correlations would help in predicting flow velocities, 
pressure drops and heat transfer coefficients. 


New experimental techniques are also called for 


especially in the measurement of low-velocity non-is ©thermal 
flow, the direct measurement of heat transfer coefficients, 
mixed convection flows, measurement of semiconductor device 
temperatures and flow measurements by infrared and optical 
techniques. 


Literature Survey 

One of the earliest books on cooling of electronic 
equipment is by Steinberg (1980) . The book deals with 
empirical rules that must be observed in designing instru- 
ments. However it does not have thermal analysis that would 
predict regions of temperature peaks called as hot spots. 
Much of the experience in this area has been developed in 
the industry and this is available in the literature in the 
form of books (Kakac et al (1985), Kakac et al (1987); Krans 
and Bar-Cohen (1983); Aung (1983)), 

Work has been done on the problem of discrete heat 
source mounted on a wall. Numerical and experimental results 
are available on this topic. For symmetrically heated, 
isoflxax surfaces, local heat transfer coefficient have been 
measured by Siegel and Nora's and both fully developed and 
entry region solutions have been obtained (Aung, 1972; 

Aung et al 1972) • Experiments performed by wiotz and 
Stutzman (1982) yielded a correlation for the maximum plate 
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temperature which is in good agreement wrth the predictions 
of Aung et al (1972). 

Conditions corresponding to heating at one plate, with 
the other plate insulated, have received special consider- 
ation. Numerical simulations (Miyatake and Fujii, 1972; 
Miyatake et al; 1973) as well as experiments (Sparrow et al 
1984b; Miyamoto et al, 1986) have been performed for both 

isothermal and isoflux conditions at heated plate. For 

14 

modified Rayleigh numbers upto 3x10 . Miyamoto et al 

(1986) observed transition to turbulence in upper portion 
of the channel, with an accompanying enhancement of the 
local Nusselt No. at the heated isoflux surface. 

Afrid and Zebib (1989) have discussed natural con- 
vection air cooling of a heated component mounted on a 
vertical wall. They have carried out a numerical study 
natural convection air cooling of single and multiple 
uniformly heated devices. They have used a two dimensional, 
conjugate, laminar flow model. Effects of component thick- 
ness, the spacing between the components, nonpowered com- 
ponents and highly powered components have been studied. 

A qualitative flow visualization and quantitative 
experimental study has been done by Tsang-Yuan Lin and Shon- 
Shing Hsich (1989) on natural convection heat transfer in 
vertical channels with asymmetrically placed discrete heated 
ribs. Both opposing and assisting flows configurations have 



been considered. Based on the analyses of photographs and 

interf erograms, it is suggested that turbulent flow should 

be expected when the local modified Rayleigh number is in 

7 9 

the range of 1.29x10 to 9.43x10 • The heat transfer data 
and flow visualization photographs indicate that the 
present rib geometry and stratification are two major reasons 
influencing the temperature of the heated ribs. 

Ravine and Richards (1988) have studied natural 
convection from a discrete thermal source on a channel wall. 
This paper presents the results of an experimental study 
of natural convection heat transfer from a discrete thermal 
source located on the wall of a vertical channel, 

decent studies in natural convection in rectangular 
enclosures with protruding heaters on one vertical wall are 
by Kelleher et al (1987) . Jjee et al (1987) and Liu et al 
(1987). The first two papers report experimental (Part I) 
and numerical (Part II) studies of geometric effects of single 
two dimensional heaters mounted on one adiabatic vertical 
wall of an enclosure. 

Work has been done on the problems involving forced 
convection. Aung (1983) presented heat transfer meas^are- 
ments for laminar flow past a backward facing step where 
the step walls are maintained at a constant tenperature. 
Results are given for the local heat transfer distributions 
as well as the average value for all the step surfaces. 


Shou-Shing Hsieh, Huei-Janshih and Ying-J’ong Hong 
(1988) have studied laminar forced convection from surface 
mounted ribs. The study includes numerical and experimental 
aspects on low speed forced convection heat transfer near 
two dimensional transverse ribs in which the walls are 
held at a uniform temperature. 

Early experimental studies of forced convection 
heat transfer from flush mounted sources are those of 
Baker (1972, 73), who considered silicon oil and R-113 in 
parallel flow over small chiplike heaters. Single-phase 
forced convection coefficients are seen to increase signi- 
ficantly with decreasing surface area and substantially 
exceed results associated with nucleate boiling as well 
as predictions based on two-dimensional boundary layer 
theory. Studies related to heat transfer from an array of 
discrete isothermal source flush mounted on one wall of a 
rectangular channel have been performed by Incropera 
et al (1986) . 

A common electronic package is one that involves 
forced air cooling of protruding components on a printed 
circuit board. Many such boards are stacked in a direct- 
ion normal to the flow, forming parallel plate channels 
with the component covered board facing the smooth side 
of an adjoining board, f^'espite enormous variability of 
geometric conditions, systematic attempts have been made 



to determine component coefficients and to develop methodo- 
logies for determining the effect of a component location 
in an array. In an early study Duller and Kilburn (1981) 
obtained heat transfer data for a single rectangular compo- 
nent and successfully correlated the data using a hybrid 
characteristic length based on features of both the flow 
and component geometries. In subsequent experiments. Sparrow 
and co-workers considered an array of component and system- 
atically examined the effect of row number, missing and 
irregular components and barriers on the average convection 
coefficient for a single component in the array (Sparrow et al; 
1982b, 1984a; Sparrow and Ots, 1985) • Lehmann and Wirtz 
(1985) also considered the effect of rib spacing and channel 
height by placing ten unheated, repeating ribs upstream of a 
heated, two dimensional rectangular component. The effect 
of protruding ribs on flow and heat transfer in channels has 
been studied by Han et al (1984), Tanasawa et al (1984), 

Webb and Ramadhyani (l985) . Ichimiya (1987) and Kelkar 
and Patankar (1987). 

Under mixed convection conditions the flow is driven 
by an externally imposed presstare gradient as well as by 
buoyancy forces, ^f the channel is vertical, buoyancy acts# 
to augment or retard the flow# the channel is horizontal 
and heating occurs at the bottom of surface and buoyancy 
may induce a secondary flow. Heat transfer enhancement is 
most pronounced when the channel is horizontal and heat is 



transferred from the bottom surface. In experiments per- 
formed for laminar flov; between asymmetrically heated 
parallel plates (Osborne and Incropera, 1985) convection 
coefficients at the bottcxn plates were found to exceed 
those corresponding to pure forced convection. Enhance- 
ment was due to buoyancy -driven flow which replaced 
warmer parcels of fluid ascending from the plate with 
cooler fluid descending from the main flow* 

Studies of mixed convection heat transfer from 
discrete sources are limited to the work of Kennedy and 
Zebib (1983) and Tomimura and Pujii (1987). For small 
heater strips flush mounted to the walls of a horizontal 
channel, Kennedy and Zebib observed the formation of 
logitudinal and spanwise vortices in air due to heating 
at the bottom and top surfaces respectively* 

It is clear from the above siorvey that the ele- 
ctronics industry has evolved rules for thermal design 
based on experience but not analysis. On the other hand, 
analysis by research groups has been restricted to study of 
forced, mixed and free convective flow past a single or 
distributed heaters on a flat horizontal or vertical surface. 
A few experiments exist to describe flow in the vicinity of 
the wall heaters. For arrays only experiments exist, though 
these deal with only the overall thermo-hydraulic behaviour. 


The problem of several vertical surfaces adjacent to 
each other having a distribution of heaters on them has not 
been studied at the level of details needed to evolve design 
rules for electronics applications. Certain other features 
such as the effect of protrusion on the formation of hot 
spots have not been systematically investigated. There does 
not appear to be any simple model for the prediction of 
thermal behaviour of an array of circuit boards. 

Scope of ^re_se_nt. ^ Work 

Work has been done on the problem of natural conve- 
ction cooling of discrete heat source mounted on a wall. 
Little work exists on the problem of protruded heaters. 

The standard solutions available for the problem of flush 
mounted heater have been found to be inapplicable. 

The present work aims at determining the flow patterns 
and temperature field in arrays of heat generating circuit 
boards. The circuit boards are modelled experimentally 
and the relationship between maximum temperature and regime 
of flow (such as free, mixed or forced convection) are 
sought to be determined. The effect of wall protrusion on 
flow and heat transfer are studied. Optical methods are 
useful in visualizing flow and thermal field and obtain- 
ing some of the relevant parameters. They are non-intrusive 
in nature and have the advantage of calculating local values 
at a particular point. In the present work interferometry 
is used to get qualitative information about the flow field 



and to a cross check the results obtained by conventional 
techniques. In particular visualization of isotherms in 
the space between vertical flar surfaces with a distri- 
bution of discrete wall heaters using a Mach-Zhender inter- 
ferometer is one of the objectives of the present study. 

The work reported here is experimental in nature in 
an area which is vital to the electronic industry. The 
heat transfer problem is expected to become a serious 
limiting factor, once the industry starts using large scale 
integrated circuits and transputers. The work will be 
helpful in increasing the knowledge-base existing to 
address the problem of efficient heat removal from 
electronic systems. 
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Chapter 2 

EXPERIIIEMPAL SETUP AISID INSTRUMENTATION 

The problem of interest here is the determination of 
heat transfer rates in circuit boards geometries since the 
maximum temperature reached any where in the system depends 
on them. Heat transfer is strongly controlled by buoyancy- 
driven motion of air in the space enclosed by circuit boards, 
by a superimposed draft of air or a combination of the two. 

In the absence of external flow, air movement arises from a 
density differential between the light air in the heated 
area and heavier air of the ambient. The chip is mounted on 
a circuit board like a protruded heat source mounted on a 
wall and the circuit boards are stacked parallel to each 
other in a cabinet. This can be experimentally modelled by 
a heat source in the form of heated metallic strip mounted 
on an insulating sheet plates. A second passive plate is 
kept parallel to the active board to simulate the circuit 
board array, 

1 

L 1 

r 

s 

y 

Fig, 2,1 : Natural Convection Single Chip 2-D Configurations 




The parameters governing the problem are the ratios 

L/h and H/o# the heating rate of each of the strips mounted 

on the vertical plates, the relative positrons of heaters 

on adjacent boards and the extent to which they block the 

flow. In the present study measurement and analysis are 

three 

carried out for a two-dimensional heater at steady state. 

Design of Experiments 

The experiments conducted can be divided into two 

parts . 

(a) Natural Convection Experiments 

(b) Forced Convection Experiments 

(a) Natural Convection Experiments ; Experiments have been 
carried out for both single chip and multichip heaters. In 
the experiments with a single chip the heater consists of an 
aluminium strip. High quality nichrone wire is used to 
heat the aluminium strip. The resistance of nichrome wire 
is almost a constant for any change in temperature. The 
aluminium strip is mounted on a bakelite sheet having a 
groove cut on it. The strip xs attached to the sheet with 
the help of brass screws. The nichrane wire is wound on a 
thin mica strip and is kept between the aluminium strip 
and bakelite sheet in the groove. To avoid direct contact 
of high temperature nichrome wire with aluminium strip or 
bakelite sheet, additional mica strips are placed between 



nichrome heater and alximinium strips and nichrome heater and 
bakelite sheet. To prevent the heat transfer to the rear 
side of the bakelite sheet it is properly insulated on the 
back side with a thick thermocoal sheet. 


Thermocouple connections 



?ig,2,2 : Single chip two dimensional heater 

To measure the temperature of heated strips Chromel- 
Alumel thermo-couples are used. The thermo-couple 
holes of size 2 mm diameter and 3 mm deep were drilled on 
the different portions of aluminium strip and the thermo- 
couple bead is fitted into them. The bead is kept in close 
contact with the strip. 

The passive wall consists of smooth perspex sheet* 
The bakelite sheet with the heater is mounted ©n a stand 
and the passive perspex sheet is kept parallel to it. 

The distance between the two can be varied by moving the 
perspex sheet. The whole structure makes an enclosure. The 




top and bottom sides of the enclos\ire are closed with some 
insulating plate. This setup generate a problem of 


Perspex sheet 


adjustable 

Fig, 2, 3 : A closed enclosure 

conduction of heat from the heater to the bakelite sheet on 
the sides A and B of the heater. This phenomenon gives a 
distortion to the assumption of discrete heater mounted on 
a vertical wall. Moreover the problem of heat convected 
from the ends C and D of the heater violates two dimension- 
ality of the problem. To avoid the heating of bakelite 
sheet and three dimensional heat transfer effects the test 
cell is modified. This is shown in Figure 2.^ , The 

aluminium strip with nichrane heater is mounted on the 
central bakelite strip and two more bakelite strips are 
joined to the central strip with a gap of 3 ram kept between 
the strips. This air gap between the strips reduces heat 
conduction of upper and lower sides (A and B) of the 
heater. 


Joining strip_ 


balcelite sheet 
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Bakelite 

piece 


A 



B 


Bokelite 

piece 


Air gap 


FiQ 


: Prevention o2 co-:uctive heat transfer to 
supporting sheet 


To avoid the three dimensional effect, two bakelite 
pieces were fixed on the both the ends (C and D) of the 
aluminium strip. This prevents the air from coming in 
contact with the heated strip and hence minimizes heat 
transfer from the ends. The multichip heater test cell 
consists of three copper strips mounted on the bakelite 


strips. Nichrome wire is used as a heater as in the 
single chip test cell. In this setuo air gap is used 
between the two consecutive bakelite strips and bakelite 
pieces are used to avoid three dimensional effects. Chromel— 
Alumel thermocouple is used to measure the chip temper- 


ature • 


To study the effect of a neighbouring wall^, a 
perspex plate is kept parallel to active bakelite plate. 

The current through the three heaters is controlled by 
means of a variac. Each heater has independent connections 
to control the supply voltage. 




Fig. 2.5 

To study the effect baffling the flow inside the 
enclosure a thin metallic strip is used. The metallic strip 
can be mounted anywhere inside the test cell by pasting it 
on the wall. The length of metallic strip is chosen to be 

Metallic strip 


Fig. 2.6 

equal to the length of heater. Larger lengths could block 
the mean flow and smaller lengths could be ineffective. 

Forced Convection Cooling 

A test cell has been fabricated to study the problem 
of forced convection cooling of electronic circuit boards. 
The test cell fabricated was vertical in orientation so as 
to have combined forced and mixed convection effects. It 
can be classified as an open circuit, and low speed wind 
tunnel. Figure 2.7 shows the overall view of the test 





cell. The material used to construct the v;ind tunnel is 
perspex. The main parts of test cell are: 

a. Entrance Cone 

b. Honey comb and filter 

c. Test section 

cj. Contracting adapter 

e. Centrifugal blower 

f. Rotameter. 

a. Entrance Cone ? Entrance cone nas two parts, settling chamber 
and contraction cone (Fig. 2,8 ), Flow from the filter 
enters the settling chamber and becomes stable before enter- 
ing the contraction cone. Small eddies formed by the screen 
and honey comb gradually decay. 

Contraction Cone s Plow from the settling chamber enters the 
contraction cone. The contraction cone is a gradually conver- 
ging cone. The decrease in cross-section takes place very 
smoothly. It decreases the free stream turbulence, increases 
the velocity and keeps boundary-layers thin. The contraction 
ratio in this test cell is 4. 

b. Honey Comb ? Honey combs are used to suppress the level of 
free stream turbulence that cones as input flow. The 
transverse component of velocity is suppressed when the flow 
passes through honey comb. Due to growth of shear instability 
and subsequent randomization of eddies, turbulence is 



generated at the tips of honey comb cells. But this high 
frequency turbulence decays rapidly with the distance result- 
ing in a net suppression of free stream turbulence by honey 
combs. The flow observes a small pressure drop while passing 
through the honey comb. The honey comb used in this test 
cell has a dimension of 500 mm x 250 mm x 75 mm. 

Filter s A fine mesh screen is placed at the downstream of 
honey comb. It helps in the suppression of turbulence 
level. 

c* Test Section ; After the contraction cone the flow enters 
the test section. Sane length is provided for the turbul- 
ence to decay. The area of the test section is 250mmxl25mra. 
It has openings for locating static taps and inserting other 
probes. The test object that is single chip or multichip 
heater can be placed inside the tunnel. On one wall of 
test cell a recess is made to accommodate the heater which 
is mounted on the bakelite sheet. On the back side of the 
heater an air gap is used to prevent energy loss to the 
ambient. The heater assembly can be mounted on the wall of 
the wind tunnel with the help of screws. The connecting 
wires for heater and thermocouple are taken out from the 
back side of heater through a hole mode in the perspex 
sheet • 

Inside the test section a perspex plate is placed 
parallel to the heater assembly to study the effect of the 
neighbouring wall. The distance between the heater assembly 



and this perspex plate is kept adjustable. The end of this 
perspex plate facing the incoming flow is wedge-shaped so 
that the mean flow is disturbed to a minimum. 

On the sides of the test section slots are made to 
have a provision for optical windows. To do the interfero- 
metric studies optical windows are required. Ordinary glass 
can not be used as windows for interferometric studies 
because it produces its own effects on light beam. A high 
quality polished glass is used as an optical window. The 
size of the glass is 50 mm x 50 mm. It is fixed in a 
perspex frame and can be slided over the slot. When inter- 

ferometry is not done, the slots can be covered with an 
ordinary perspex sheet. Figure ( 2,9 ) shows the test 

section of wind tunnel. 

Contracting Adapter ; 

It joins the test section with the blower opening. 
Since the dimensions of test section are quite large a 
contracting adapter is made. The area contraction ratio for 
adapter is 15 and height is 40cm.At the top end of the adapte 
piece of metallic pipe (diameter 2.5 cm) is fixed. This pipe 
is connected to the blower by a PVC flexible pipe, 

. Centrifugal Blower s The blower sustains the velocity up to 


0.65 m/s in the wind tunnel. The specifications are given 



below* 


Specifications 


Outlet 

Voltage rating 
Current 
Dischar ge 
Power rating 


- Radial 

- 240 V 

- 3.5 A 

3 

- 20 m /min 

- 750 W 


f . Rotameter ; Both pitot tube and rotameter are used to measure 
the velocity in the wind tunnel. Since the cross-sectional 
area of wind tunnel is large the maximum velocity reached 
with above mentioned blower was 0.65 m/s. In electronic 

oj 

equipment the air velocities encountered are^ similar order. 
Measurement of velocity as low as 0.15 m/s with a pitot tube 
is difficult. To make the measurements possible at such low 
velocities a low discharge rotameter is used. The rotameter 
capacity is 10 litres/sec. which corresponds to approximately 
0.5 m/s velocity. Hence it was possible to measure the 
velocities on the lower side with the help of the rotameter. 
When measuring higher velocities the rotameter is kept out 
of circuit. 


Instrumentation 

The experiments involved instruments for the measure- 
ment of current, resistance, voltage, temperature, pressure, 
velocity, turbulence. For interferometry a Mach-Zhender 
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interferometer having laser beam as a source of light is used. 

Measurement of Electrical Quantities 

Measurement of current, resistance and voltage was 
carried out with digital multimeter. The multimeter used 
were: 

(a) Meco model 9A digital multimeter 
Specifications : 

Operating Principle : Dual Slope integration 
Range : 100 V DC, 750 V AC and lOA AC/13C 

Power Supply : Two polarity 

Display : LCD 

Accuracy : AC Voltage: + 1*5% of input 

+ 8 counts 

DC Voltage: + 1% of input 

+ 4 counts 

Resistance: + 0,7% input 

+ 4 counts 

Current : + 2% of input 

"” + 4 counts 

A more sensitive and versatile multimeter was tised to 
measure the millivolts generated by thermocouples in some 
experiments. 

Make : Hewlett Packard 
Model : 3457 A, 

1 1 

It has seven functions with 3 ^ to 6 ^ digits of resolution 
extendable to 7 ^ digits at reading rates from 1 reading 
every 2 sec, to 1350 readings/sec. It has seven functions - dc' 
acV, acl, del. Ohms, frequency and period. 


Operating Temperature: 


J4 

0 to 55°C 

Accuracy : DC voltage - 0.005054 -t 7 count 

AC voltage - 0.254 + 280 counts 

DC current - 0.0754 + 604 counts 

AC current - 0.3554 + 280 counts 

Resistance - 0.0045 + 34 counts 

Frequency - 0.01% 

Variacs are used to step down the voltage. To stab- 
ilize the supply, two variacs in series are used. Voltage 
stabilizers of the servo-type are used to prevent under and 
over the voltages in the input. To protect the sensitive 
instruments from overloading a ferromagnetic type voltage 
stabilizer is used. 

Measurement of Temperature 

Thermocouples are used for the measurement of temper- 
ature. The thermocouples used are 18 gage K type thermo- 
couple made of Chromel and Alumel. The thermocouples are 
formed by gas welding. Each thermocouple is individually 
tested against an standard thermocouple having ice point 
as a reference junction. Variation upto +0.1% are possible 
among the thermocouples over the range of temperature covered 
in this study. The thermocouples were connected to the 
30 channel hybrid recorder. The description of which is 



given below 


Hybrid Recorder 

It has a capacity of raking input data from 30 cha- 
nnels simultaneously. It can be used for various types of 
thermocouples over a wide range of temperature. The reading 
can be displayed directly in the terms of emf generated or 
can be displayed in terms of temperature. It uses internal 
compensation instead of the ice point as a reference junction. 
There is no need to have ice point as a reference junction 
or adding the mV corresponding to room temperature eveiry 
time. It can display the variation in temperature with time 
in the form of a graph. The data can be stored in the form of 
a graph. The data can be stored in the memory and can be 
printed. The instrument can be connected to a computer over 
a GPIB interface. 


Specifications 
Thermocouple type - K 

Measurement range 

jf^esolution 

Accuracy 


- 200 to 1370°C 
0.1°C 

+ 0.07% reading + 0.8°C 


Thermocouple type - 
Measurement 
Resolution 
Accuracy 

Ambient temperature 


J 

range =* -2 00 to 12 00^^0 

= -0.1°C 

= + 0.07% + 1.0°C 

23+5*^C, Humidity 80% or less 


36 

. o ,o 

Measurement temperature unxt - C/ 5’ selectable 
Reference contact correction - Internal/E'xternal selection 

possible . 


The other temperature measuring device used was digital 
thermometer. 

Make and Model - Therm 2280-3 
Temperature range- 


Ni(Cr) - Ni(K) 
Fe-Cu-Ni (U) 
PtRh 1G-Pt(8) 
NTC 

Resolution 


-200 to ISVO'^C 

-200 to 600 °C 

0 to 1760°C 

-30 to 80°C 

0,1 K (thermocouples) 


It is easily portable and gives the readings very accurately. 


Pressure and Velocity Measurement 

For pressure and velocity measurement a pitot tube 
connected to a sensitive manometer was used. This mananeter 
can measure very low velocities and hence flow measurements 
could be performed in the mixed convection region. 

Make - Fur Controls Limited 

Maximum AP - 19.99 pascals 

Velocities below 0.2 ra/s are difficult to measure with a pitot 
tube. To measure the velocity in this range a rotameter is 
used. The discharge of rotameter varies from 0 to 10 lit/sec. 
The maximum discharge corresponds to a velocity of 0.5 m/s in 
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the tunnel* 


Turbulence Measurement 

While testing the wind tunnel the turbulence level 
was measured with the help of an rms voltmeter. The rms 
voltmeter gives the voltage corresponding to fluctuating 
component of velocity. The dc voltage corresponding to 
mean velocity can be taken out from manometer. The ratio 
of rms voltage to the mean voltage gives the turbulence 
level. 

Specifications of rms voltmeter used are -- 

Make - Hewlett Packard 


Frequency 

response 

Range 


- 10 Hz to 10 MHz 

- 1 mV to 300 volt (l2 ranges) and 
-72 to +52 dBm. 


Interferometry 


The interferometric study was carried out using a 
Mach-Zhender interferometer. The Mach-Zhender interferometer 
is a popular instrument used in experimental studies of heat 


transfer in fluids 
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Fig, 2.10 

Laser is used as a coherent monochromatic light source. 

The laser beam passes through a spatial filter. Spatial 
filter expands the beam which is then passed through a 
convex lens which makes it parallel. The beam goes through 
the beam splitter. The beam splitter is semipolished 
mirror which allows 50% of tight to go through it and 
reflects remaining to a mirror. Unreflected part of 
beam passes through the heated test section and then falls 
on mirror which reflects it at 90°C. The other half of the 
beam reflected from beam splitter also falls on a mirror 
which reflects it at 90^, Finally the reference beam and 
test beam combine at the second beam splitter where the 
interference patterns are formed. 
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These interference patterns can be recorded by a 
field camera and the analysis is done. Quantitative experi- 
ments are possible when the temperature field is uniform 
parallel to the beam direction. Otherwise the test beam 
averages the temperature field as it traverses the test 
cell and only qualitative information can be obtained. 

The fringe pattern obtained can be directly scanned by 
a ccd camera coupled i/ith an image processing equipment. 

The film can be recorded in the memory of computer and 
can be obtained on the screen of monitor. The analysis 
can be directly done through software. In the present 
study all the fringe analysis is done manually using a 
low magniLication microscope. 

Testing 

Testing in natural convection experiments involved 
testing of thermocouples and testing of the heater element. 
Each thermocouple is cleaned and tested against the ideal 
thermocouple with an ice junction as reference point. A 
thermocouple showing a widely different response is 
rejected. Before experiments the heater element is tested 
several times. It is ensured that the bakelite plate 
withstands the temperatures encountered during the experi- 
ments. This can be as high as 150°C. Proper insulation is 
provided on the rear side and it was made sure that no 
significant rise in temperature takes place at the surface 
of insulating material exposed to the ambient. 



In forced convection experiments fche testing of the 
flow apparatus is important. The test cell is examined for 
the leakages, pressure drop, turbulence level and velocity 
profile. The openings, slits or holes, if any were sealed 
properly with the foam packing and tape to prevent any leakage 
of air from outside. 

When tested for pressure drop it is found that the 
pressure drop was very low across the honey comb. Most of 
the resistance to flow occurred in the converging adapter 
and the pipe connecting the adapter to the blawer. 

The turbulence level was measured using pitot tube 
connected to a manometer and a rms voltmeter. The maximum 
turbulence level was 2.5% at the low velocity of 0,24 m/s but 
as the velocity was increased the turbulence level was found 
to decrease. At the high velocity of 0.67 m/s the turbulence 
level was 1.0%, These values of turbulence levels are accep- 
table at such low velocity ranges. This is because for weak 
flows it is difficult to protect them from the external 
influences. This turbulence level includes the flow distur- 
bances due to the movable plate, fluctuations in voltage 
supply and the vibrations caused by the blower. The 
table gives the value of turbulence level with velocity. 

The turbulence level can be seen to be constant in the 
direction perpendicular to wall* 
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As tihe velocity changes the boundary layer thickness changes. 

At loiter velocities the boundary layer thickness is large and 
the heater is under the influence of boundary layer. It is 
difficult to estimate the effect of this var-^ing velocity pro- 
file on heat transfer. Heat transfer rates are also different 
over the different parts of heater. Here v/e have not taken 
the velocity profile and boundary layer effects into consider- 
at ion and it is fairly accurate to assume the velocity in 
the centre of channel as the reference velocity for calculating 
the i^eynolds number over the heater. Tne reason for this is 
explained belovr. 

The front face of neater is the region where most of 
the heat transfer takes place. The upstream region of heater 
is, to some extent, a stagnant region not affecting the overall 
heat transfer much. Enveloping of this region inside the 
Doundary layer does reduce the heat transfer but a corres- 
ponding increase takes place on the front surface of heater 
because of increase in velocity as the area reduces in that 
region. The affect of upstream face of heater being subjected 
to a reduced velocity is thus compensated. Moreover oecause 
of protrusion the front side of heater is subjected to the 
velocities prevailing in the central portion of channel. Due 
to the combined effect of these two factors it is reasonable 
to assume that the velocity in the central part of the channel 
IS the velocity having maximum influence on beat transfer. 
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This is also the reason of choosing length L oi front side 
as the characteristic length in this ^j'ork. 

The shape of the velocity profile does not seem to 
affect the heat transfer raucn. This has been shown by Hsieh 
et al (1988) in their numerical work, '^e have also seen 
similar effects in the numerical calculation of the Musselt 
over a flush mounted heater. These results have been presen- 
ted in the Chapter 4, 


Turbulence Measurement 


Probe i^ear the Wal 1 


Velociry 

m/s 

bean 

DC voltage 

rms voltage 

% Turbulence 

0,24 

0.23 

0.006 

2.60 

0.33 

0.30 

0.007 

2.30 

0.37 

0,34 

0.007 

2.00 

0.42 

0.38 

0.0075 

1.97 

0.47 

0,43 

0,008 

1.86 

0,52 

0.45 

0.0065 

1.41 

0.57 

0,50 

0,007 

1.40 

0.61 

0.53 

0.006 

1.13 
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Probe in the Middle of the Test Section 


Velocity 

m/s 

Mean 

DC voltage 

rms 

voltage 

% Turbulence 

0.28 

0.24 

0.006 

2.5 

0.32 

0.28 

0.008 

2.85 

0.40 

0.34 

0.006 

1.76 

0.47 

0.41 

0.0065 

1.58 

0.51 

0.45 

0.007 

1.55 

0.55 

0.50 

0.006 

1.2 

0.60 

0.54 

0.0065 

1.2 

0.64 

0.56 

0.0055 

0.98 

0.67 

0.59 

0,006 

1.0 


The maximum temperature encountered in the experiments is 
nearly 150°c. So a wide range of has been covered. 

In natural convection experiments with a heater of AR = 3# 
the Rayleigh number varies from 2000 to 70000. With 
smaller of aspect ratio this range is narrowed. The 
maximum value of input heat flux is limited by the maximum 
temperatures reached and the power drawn. 

In forced convection experiments with larger heater 
the range of Reynolds number covered is from 0 to 1050. With 
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smaller heater the range is reduced from 0 to 600. Maximum velo- 
city encountered is 0,65 m/s. With the protrusion mounted on 
the opposite wall the resistance to flow is increased and the 
maximum velocity reached reduces to 0.5 ra/s , 

The fluid Prandtl number varies from 0.7 to 0.3, Two 
cases for aspect ratio 2 and 3 have been considered. 




: Test cell f< 
experiments 


Figure 2 U 



forced convection 



Fig 2il Single chip natural 
path of test beam 


ce/f 
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Fig. : Fringe pattern obtained on screen 
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Chapter 3 

Ponr’V.n TTRC AI?D t iATA REPUCTl^ 


we describe belwj the procedure for natural conve- 
ction experiments, procedure for forced convection experi- 
ments and data reduction for interferometry. 


Natural Convention Experinient^ 

T 'hfahind the aluminium chip is 

The heater located Penxna 

supplied voltage through a variac connected to a servo 

stabilizer. The supply voltage to the heater is bept 

stable to within +»• The heater is allowed to reach the 

steady state. The steady state is reached in approximately 

o u rr «i-;.rt steady state is assumed to be 

3 hours from a cold start* » s 

. .^r^cr'a•^-nr•p o£ the chip is less 

reached when the rise m temper 

than 0.5% over a half-an-hc»r period. To accommodate the 
yariatiohs in input supply and variation in response of the 

individual thermocouple temperatures and voltag 

» X.U 5 rcixn'o.ts intoxvals# 

averaged at steady state with A 

• 4- • snd tQinpsr-" 

This averages out minor variations, power inp 

h’ noi-riimented with six thermo- 
ature with time. The strip is instrumentea w 

couples. The temperature of the strip is taken as the 

average of all the thermocouple readings. This compens 

^ 4 - 1 -.Q einrface o£ chip* It can 
variation in temperature over t 

^ •^-r^ -hP at a uniform mean temperature* 
be subsequently assumed to be a 

1- v.c4-iTro is recorded with the 

At every step the room temper 

.n. . ^ >11 the experiments are conducted in 

digital thermometer* All tne exp 



an airconditioned room maintaining a fairly constant room 
temperature. The test cell is protected suitably from the 
external air currents* 

The procedure for multichip experiments is similar to 
the single chip experiment except that in the case of multi- 
chip heater the temperature of all the chips is recorded 
and the power supplied to the chips is controlled indivi- 
dually* 

Forced Convection Experiments 

In forced convection experiment the heater plate is 

mounted in the flow apparatus in a vertical configuration • 

The electrical power is once again supplied to the heater 

through the servo-stabilizer variac. The heater is then 

allowed to reach steady state. All the readings at steady 
state corresponding to free convection are recoded. The blower 

is then switched on inducing a flow in the wind tunnel* 

The blower is made to run at constant speed by controlling 
its input voltage. This corresponds to a constant mean 
velocity in the test cell. The temperature of the chip 
drops due to an increase in heat transfer at its surface* 

The chip is again allowed to reach steady state. This 
takes about 2 hours . The temperatures are recorded at 
steady state as in the natural convection case by averaging 
over a period of time. Experiments are conducted for 
various values of the input heat flux for a range of air 
velocities* 



Interferometry 


The interf ercroeter is first adjusted in the infinite 
fringe mode or the wedge fringe mode depending on the 
requirement. The test enclosure mounted on a stand is 
placed in the path of the test beam. As the strip gets 
heated a fringe pattern is obtained on the screen. The 
photographs are taken during the transient evolution of 
free convection. The photographic film can be developed in 
the laboratory itself. The fringe pattern is recorded on 
an image processor as well. The image can be enlarged and 
viewed. The fringe pattern has been analyzed in this study 
under a low magnification microscope. The microscope has 
a traversing mechanism in two perpendicular directions. 

The distance between two fringes or the slope of the fringe 
can be measured and the heat flux can be obtained. Nusselt 
number obtained by interfercxnetry can act as a cross check 
for the Nusselt number obtained by the energy balance method. 
The behaviour of flow and temperature field is also studied 
in the vicinity of the heated strips using the fringe patterns. 

Data Reduction 

The resistance of every heater is measured before 
experiment. Heat input is given by 



where V is mean voltage applied (volts), R is resistance 

2 

of heater (ohms). Let A (m ) be the area of strip exposed 
to ambient. The heat flux q is 
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q = vf W/ra^ 

AR 

Let the temperature of heated strip be and the temperature 
of ambient be foo* 

Heat transfer coefficient h is determined as. 


h = W/m^K 

Let L be the characteristic length of heater. It is the 
strip length vertically in the present study. 

Let k be thermal conductivity of air. The Nusselt 
number is calculated as 


The Rayleigh number is calculated as 

3 


Ra 




va 


c 




/ fj w ^ 


1 

^ ^ tQ 

113 ^.' 


Here p is the density of air 


AT = 

P = 


'C OO 

l^^+273) 



y = Kinematic viscosity of air (m /sec) • 

^ = is the thermal diffusivity of air (m^/sec). 

(All the fluid properties are calculated at the mean tem- 
perature + "^room 5’urther g is the gravitational 

2 

constant (= 9.81 m/sec ),and L is the characteristic length 
of the heater. 


In forced convection experiments Reynolds number is 
calculated along with Nusselt number and Rayleigh number. 

A calibration curve is drawn between the velocity in the wind 



tunnel and the voltage applied to blor/jer. The velocity in 
the v/ind tunnel is measured \,Tith pitot tube connected to 
the manometer. Once the calibration is done the pitot tube 
IS taken out. 


Reynolds number is calculated as 


V 

iSJusselt number and Rayleigh number are calculated as in 
the previous experiments. 


Calculation o£ Nusselt Number by Interferometry 

When the initial setting of interferometer is in 
infinite fringe mode then each fringe obtained on the 
screen represents an isotherm. The change in temperature 
for each fringe width spacing is given by 



Yt is refractive index of medium, \ is the wavelength of 
light source (He-Ne laser source in this case) . 1 is 

length of test section parallel to light-beam. In air 
= 0.927 X IQ-V^C. 

For a He-Ne laser X = 632,8 nm. It can be shown 

that the temperature per fringe shift is 

T = 0.67/1 °C/ra. 

If e be the fringe spacing then the local heat flux is 
given by 




Nu 


0.67 

1 (t - t ) 
' c oo 



This method is suitable when the number of fringes are 
small i»e« temperature gradients existing are not very 
large. When the temperature gradient is large near the 
hot surface then many fringes are formed at the surface 
of heated strip. For the calculation of heat flux at 
the surface the fringe spacing closest to wall is meas- 
ured. In the case of high temperature gradient the distance 
is too small to be measured or in other words it becomes 
difficult to distinguish among fringes. To overcome 
this problem another approach is used. 



54 




' i 9 » "3*1 

The temperature at the hot surface rs known (T^)* Seme 
fringes away from the heated wall are chosen. The fringes 
are chosen at a point where the fringe spacing can be 
measured conveniently. Since each fringe width corresponds 
to a fixed temperature difference, the temperature corres- 
ponding to fringe A, B and C can be found out by counting 
the number of fringes formed from the surface of strip. The 
temperature T^, Tg and T^ are known. The distance of the 
fringes A, B and C is measured from the wall* T^ be the 
temperature at a distance Tg be the temperature at a 
distance y 2 # T^ be the temperature at a distance 'y 3 . 


A 

curve < 

2an 

be fitted 

i passing 

through 

these points 

Let T = 

^h 

ay 

■ + by^ 

3 

+ cy 

then 



= ^h 


a Yl + 

by2 + 

cy^ 

(1) 


= ^h 


a y 2 + 

i>Y2 + 

3 

(2) 

Tq 

= ^h 


a Yg + 

^ 2 

bY3 + 

3 

cyg 

(3) 



This is a system of simultaneous equation having 
a# b and c as unknown quantities. These equations can be 
solved to ' the values of a, b and c. The heat flux 
at the surface is proportional to the temperature gradient 
at the surface and is given as. 


q 


surface 


al 

ay Y= 


■0 


= ka 


This gives the heat flxix locally at any point. To obtain 
average heat flux 


q = 

q = 



JS- 

L 


L 

/ k a dx 
o 

L 

/ a dx 
o 


Nusselt number Nu 


Nu = 


o 


oo 


dx 


‘Th - 



The value of 'a' can be integrated over the whole width 
L and average value of Nu can be found out. 
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The table 3 i presents the values of voltages and 
temperature taken at steady state. Each voltage reading 
was taken after 5 minutes. 


Mean voltage V = 10.84 volts 

o 

Mean chip temperatiire T^, = 43.15 c 

o 

Mean room temperatxire T^ = 29.67 ^ 

Resistance of heater = 69,2 Ohm 

"^c ^oo _ 43.15 + 2 9.67 

2 2 

= 1»70 

R 69.1 



36.41 


q , A = 7. 985X10’’^ m^ 

q = lil = 212.89 W/m^ 

7.985x10“ 

h = « — — = 212,89/13.48 = 15.79 

(T - T ) 

^ c oo 


Nu = ~ (k is thermal conductivity of air at 

k 

mean temperature. The value can be 
taken from table) 

k = 27o00xl0“^ W/iriK 


Nu 


Ra 

P 


14.618 


pgAT' 


va 


1.139 kg/m^ 
/ 2 

9.81 m/sec 


g 



Ch.No. V=10.85 V=10.8 V=10.81 V=10.79 V=10o8 V=10,94 V=10.95 Mean 
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(36«4l + 273) 


— 5 2 

= 1,661x10 m /sec 

“>•5 2 

a = 2.35x10 m /sec 

L = 2,5 cm 

4 

Ra = 1,69x10 


= 3o231x10”^/k 


Forced Convection Experiment 

The data for the sample calculation is given in 
the table 3 2. 

2 2 

q = W/m V(volts) is voltage applied to 

heater, V = 2 9,5 

A(m ) IS the area of strip exposed 
to atmosphere 

-3 2 

A = 7.74 X 10 m 

R = 110,8 Ohm 

L = 2 , 5 cm 

q = 121.^5) = 1014,7 W/m^ 

110.8 X 7.74 X 10“ 


Velocity = 0.15 m/s 


T, 


58.25°C 


T 


CX5 


mean 




19,1°C 

58.25 + 19,; 
2 


331.25 


A 


38.67°C 



All the fluid properties are evaluated at mean temperature. 


Nu 


hL 

k 


h = 


cr 

AT 


h = 


1014,7 


(58o25-19.1) 


= 25.91 W/m K 


Hu 


23 .87 


R,s 


pg ZiT L' 


= 45377.4 


va 


Ijji — 0,15 X 0.025 „ 225 S 

y 1.661x10"”^ 
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Calculation of Nusselt Number by Interferometry 
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we do a sample calculation for obtaining the temperature 

profile and thereby the Nusselt number through cxirve fitting. 

/ 2 

The power input q = 564.72 W/m 

Temperature of heated surface T^ = 48*7°C 

o 

Ambient temperature = 18.1 C 

Three fringes away from the heater are selected where the 

fringe spacing can be conveniently measured. The temperature 

o 

change per fringe shift is 3.2 c. 

Temperature corresponding to fringe 1 = 39.1, fringe 
2 = 35.9, fringe 3 = 32.7. The fringe width is calculated at 
6 locations along the vertical length of the heater. 

oc- 



The fringe width at these points is given below; 


Location 

mm 

X 

^1 

^2 

^3 

a 

X 

Hu^ 


mm 

mm 

mm 

mm 



1 

0 

0.15 

0.75 

2.85 

-78.76 

64.34 

2 

5 

0.375 

1.65 

3.75 

-30.04 

24.54 

3 

10 

0.45 

1.80 

3.90 

-25.53 

20.85 

4 

15 

0.525 

1.80 

3.90 

-22.45 

18.40 

5 

20 

0.525 

1.725 

3.75 

-22.59 

18.44 

6 

25 

0.30 

1.50 

3.30 

-36.82 

30.08 
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Ylf Y2 y3 fringe widths for the fringes 1,2 

and 3 respectively. 

2 3 

we have T = + ay + by + cy 

Substituting the values for T2# 'r3 and corresponding fringe 
V7idths in the above equation we get for location 1 . 


33.1 = 

48.7 


a 

(0.15) 

-h 

b 

, n2 

(0.15) 

+ 

c(0.15)' 

35.9 = 

48.7 

+ 

a 

0.75) 


b 

(0.75)^ 

4- 

a 

c(0.75) 

32.7 = 

48.7 


a 

(2.85) 

4* 

b 

(2.85)^ 

4* 

c(2.85)^ 


These three equations can be solved to give the valu« of a, 
b, and c. Similarly rhe value is obtained for all the points 
along the lengtn of heater surface. 

As we have said earlier 


M.\ 

■2>y|y=0 


a. 


The value of temperature gradient at chip surface is known 


q 


X 


■k 


'bY y=o 


Nu^ - 



L 



T 


00 


) 


— k a 

X 

(k = 27 x 10 “^ w/mk) 


The average Nusselt number is obtained by integrating these 
values over the entire length of the heater 


Nu = 





dx 


Using trapezoidal rule we evaluate the average value of 
heat flux and Nusselt nrmaber as 
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Chapter 4 

RESULTS AIID PIS CUSS I OH 


The results are presented below for the follOT,jing 
configurations ♦ 

Single chip natural convection* 

Three-chip natural convection. 

Single chip forced convection. 

Three chip forced convection. 

Qualitative and Quantitative analysis by inter- 
ferometry 

Single Chip Natural Convection 

In these experiments the variation of Husselt number 
with Rayleigh number has been studied. The chips are made 
of aluminium. In view of the high thermal conductivity of 
aluminium these can be assumed to be at a uniform temper- 
ature over their surfaces. Direct calculations show a neg- 
ligible temperature drop from the thermocouple bead position 
to the chip surface. The effect of placing of neighbouring 
passive wall has also been determined. The role of heat 
conduction to the supporting bakelite plate and the effect 
of radiation losses has been investigated and eliminated to 
a large extent. Hence the Nusselt number reported here 
corresponds entirely to the convection mechanism in the fluid 
phase* 
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The Nusselt number shows an increase v;ith the increase 
xn the chip temperature and hence the Rayleigh number. 

Table (4,1) gives data for an enclosure that is 

kept open at the top and bottom. The distance between the 
active plate and neighbouring passive plate is 3H where 
H is the height of protruded heater. Figure 4,3 compares 
the data for an enclosure that is closed at the top and 
bottom with an enclosure that is ooen on these sides- 






□ 


Figure 4,1 '~'dosec( 

Nu increases with Ra in. general but in case 

of a closed enclosure the Wusselt number obtained is less 
than that in an open enclosure. The possible ncchanisms that lead 
to this result are as follows. When the enclosure is open air 
movement is free and uninterrupted. Owing to buoyancy 
air rises from the lower side and leaves the top of the 
enclosure in a vertically upward direction* 


lit 




fit 


Fig. 4.2 : Buoyant flow through open endosurt 




i 


when the enclosure is closed at the top and bottom faces, 
the cold air from the lower portion of enclosure rises but 
can not exit from the top face. Hot air gets accumulated 
in upper portion of the enclosure and results in an incre- 
ase in overall temperature of the enclosure. The possible 
air flow paths are shown below. This is a case of combined 
horizontal and natural convection. Air is taken in from 
the lower sides of the enclosure and leaves the enclosure 
on a coiribined horizontal and vertical path. 

The absence of accumulation of heated air in the 



Fig. 4.4 : Plow mechanism in a closed enclosure 

upper portion of enclosure and a shorter, more vigorous air 
flow mechanism results in heat transfer being better in 
the case of an open enclosure. This emphasizes the need of 
proper ventilation in cabinets having many circuit boards. 
The systems having natural convection cooling should have 
sufficient space below and above the heated elements. 





Table 4,2 presents data for the closed enclosure with 
the passive plate hept at a distance of 2H from active 
plate. On conparing Tables 4,2 and 4.3, it ca?-i be seen 
that the position of passive plate hardly aff-cts heat 
transfer. The distance betti/een the two plates is not a 
critical parameter unless it is made excessively small. Ihe 
observation can be explained as follov/s. In natural conve- 
ction there is a formation of a boundary layer over the 
heater. The length of the heater parallel to this flow is 
quits small. Hence the thickness of the boundary layer is 
also small. Since all the terperature drop takes place 
within the boundary layer and the fluid outside the boundary 
layer is at the ambient temperature, the neighbouring plate 
does not disturb the boundary layer. As a result the heat 
transfer rates are not affecbed. If the neighbouring plate is 
kept very close to the heater xt decreases the KFusselt number. 
This is due to the unavailability of a sufficient amount of 
cold air in the vicinity of the chip. 

The effect of a sharp protrusion kept close to the 
heater to divert the flow towards the heater has also been 
studied. The effect of the protrusion is studied for two 
values of heat flux. The Table 4,4 shows the effect of 
protrusion placed at various positions, The different posi- 
tions are shown in Figure 4a, 
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It can be seen from rhe data presented in this table (4,4) 
■chat Nusselt number is lo'uered when the disturbing strip 
is placed in the path of the flow. The protrusion tends 
to block the flow resulting in an increased resistance and 
a reduction in average air velocity. For a given heat flux 
the Nusselt number obtained without protrusion is higher 
and hence the maximum temperature reached on the chip is 
lowered. 

Effect of Radiation from the Chip 

To estimate the amount of radiation losses the experi- 
ments were separately carried out in a white enclosure and 
a black enclosure. In a white enclosure all the radiation 
falling on the walls is assumed to be reflected back to the 
heated chip while in case of black enclosure all the radiant 
energy falling on the walls is assumed to be absorbed comple- 
tely. 

It can be seen fron the Fig. 4,5 that for given Ra, 

Nu is higher in case of a black enclosure. Here the heated 
metallic strip loses a part of its energy by radiation. This 
results in a lower surface temperature. In case of a white 
enclosure the radiation mechanism is not available and all 
the energy released goes to the ambient air by convection 
only. This leads to a higher chip temperature and subse- 
quently lower Nusselt number.The amount of losses due to 
radiation for the test cell used in this study are approxi- 
mately 1(K of the total input heat flux. Appendix A.3 




Ra 

Fig •4. 5 ; Effect of Radiation 
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the chip 



gives the detailed method of calculating the radiation loss- 


es . 

Tne Table 4.7 gives the values of radiation loss in 
a black enclosure. If experiments are done in a black 
enclosure then the data free of radiation effecbs can 
be obtained as of the measured ’Junselt number. In 

present study we have used v^hite enclosure for all the 
natural convection experiments and this correction is not 
necessary. 

Uall Conduction Effects 

As the temperature of the strip increases the temper- 
ature of the wall on v;hich the strip is mounted, also 
increases because of conduction. This changes the temper- 
ature profile over the bakelite sheet and precludes the 
assumption of discrete heater. This change in temper- 
ature profile affects the overall heat transfer rate. Air 



Fig.4.6 


rising upwards comes in contact with the wall first and gets 
pre-heated before it reaches the heater. The Nusselt number 
calculated on the basis of the assumption, that all the heat 
supplied is lost to ambient through the strip only, is 




likely to be higher than the Nusselt number in which end- 
conduct ion is present. 


To minimize the effect of conduction, experiments are 
carried out in a test cell having an air gap between the 
heated strip and the supporting bakelite sheet. Two bake- 
lite pieces are fixed at the ends of strip to avoid three 
dimensional effects. 

Air gap 



A 



Bakelite piece 


Pig, 4, 7 : Air gap to prevent conduction 


Figure 4.8 shows that the effective Husselt number 
obtained in the setup with a slit is less than that obta- 
ined v^ithout a slit. The two values differ by 5-10%, This 
result is ej^^ected as explained earlier. At Higher Rayleigh 
number this difference increases. The temperature of heater 
increases linearly with an increase in the heat flux. 

Ravine and Richards [^1988} present experimental data 
for a flush heater mounted on a vertical wall. The heater 
length is 1 and channel separation is b in their study. 






I 

d 


Pig. 4, 9 

X xs measured from the leading edge of the heater. 

For b/i > 1 they report 

Nu 

0-4 

In the present study this ratio varies from 0,9 to 
1,6. The protruded heater does not show strict similarity 
as in the case of a flush heater. The protruded heater 
generates a non-parallel flow field and hence a higher 
Nusselt nximber. 

Aspect Ratio Effects 

The aspect ratio has been defined as L/h in this work. 
Two cases for aspect ratio 2 and 3 have been considered. 
Figure 4.10 shows a plot of temperature versus heat flux 
for heaters with aspect ratios 2 and 3. Results for a three 
dimensional heater are also shown here. The heater with 
larger aspect ratio shows a higher temperature rise than for 
a heater with smaller aspect ratio. Figure 4.11 shows 

3 

variation of a Nu/(AR) with Ra/(AR) • The axes have been 

3 

taken as Nu/(AR) and Ra/(AR) so as to have the same char- 
acteristic length for the two heaters. 





Fig, 4,10 : Terrperature rise with heat flux 
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Pig* 4*11 : Aspect ratio effect (Natural Convection) 



It is evident fron the two figures that the heat trans- 
fer rate is higher in case of heater with smaller aspect 
ratio. This is because the heacer with smaller aspect ratio 
is closer to a point soiirce. The point source and a 
uniformly heated flat plate are the maximum and minimum 
limits for heat transfer. 

As the size of source becomes smaller the flow field 
becomes three dimensional, we present results for a three 
dimensional copper heater with a square cross-section mounted 
on a bakelite plate. 

Experiments show that the Nusselt number for a horizon- 
tal configuration is smaller than that for the vertical config- 
uration. Since the heater size is small the variation in 

3 4 

Rayleigh number is also small (10 ^ Ra 10 ) and flow remains 

close to laminar. The corresponding change in Nusselt 

number is also small. From Fig-* 4.10 it can be seen that 

heater 

the temperature rise for a 3-dimensional^is less than that for 
a 2-dimensional heater. Correspondingly the Nusselt number 
obtained with a three dimensional heater is considerably 
larger than for a two dimensional heater. The effect of 
placing a wall in the proximity of the heater does not signi- 
ficantly affect the heat transfer rate. 

Radiation losses are found to be insignificant in case 
of a three dimensional heater. At a heater temperature of 
53 '^C the radiation loss is found to be 1,4% of input heat 
flux. For two dimensional heater# a larger surface area and 
a large chip temperature leads to radiation losses as high as 

1 

10% of the input energy. Since the Rayleigh number for a 



small heater xs small the flow is likely to be easily effe- 
cted by external influences. The experiments dedicated to 
scatter show scatter as large as lO/a with a 3D heater. This 
e35plains the slight inconsistency in the value of Nusselt 
number with increasing Rayleigh number. A detailed scatter 
analysis is given later. 

Since the heat transfer rate for a two dimensional chip 
is the smallest, the modelling of an IC chip as a two dimen- 
sional heater will give the worst thermal performance of an 
electronic system. This will help in designing the system 
on the basis of the worst conditions. 

Multichip Natural Convection Experiments 

The aim of these experiments is to study thermal inter- 
action among the heating elements mounted on a vertical wall. 
The presence of an unpowered heater close to a powered heater 
can also affect the heat transfer characteristics. Experi- 
ments have been carried out with three similar heaters moun- 
ted on a bakelite sheet. In these experiments we have 
used heaters with a small aspect ratio to simulate closely the 
electronic conponents. The height of the protrusion 

is same as in the case of single chip heater but the width 
has been reduced to 15 mm. 







Fig. 4.12 

The three hearers are controlled separately ■» Various 
combinations of cold and heated srrips are possible. We 
have studied interaction effects for the three configur- 
ations shown below. 
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Middle Chip Heated 


In this configuration energy is supplied to the central 
chip only. Table 4.14 gives rhe value of relevant heat 
transfer parameters for this case. The Nusselt number incre- 
ases v.dth the Rayleigh number for the middle chip as shown 
in the Figure 4.14 . 

Figure 4.15 shows the variation of temperature beyond 
the ambient for all the three chips. The lower chip is found 
to be barely affected by the presence of the beared central 
chip. The temperature of upper chip increases significantly. 
This sho\TS that the presence of a powered component does 
not significantly change the temperature of a component below 
0^1 a circuit board. If the same powered component is mounted 
below, the rise in temperature of the chip above is signifi- 
cant o 

As seen in the case of a single chip heater# the open 
enclosure leads to improved cooling in the case of multichip 
heaters. The comparison between two cases is shown in 
Fig, 4,16 • Fof a given Rayleigh number the value of 

Husselt number is higher in an enclosure open fron the top. 
This is, as explained earlier, due to vigorous ventilation in 
an open enclosure. 

Lower Two Chips are Heated 

In this configuration energy is supplied to the 
chips 2 and 3. The effect of a powered conponent below and ai 
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Fig.4.l6 : Conparison 
multichip 





unpov/ersd component above has been studied. The lower chip 
is found to have a higher Nusselt nuraber than the central chi 
Phis is because for the lower chip the cold ambient air is 
available drawn from the bottom of the enclosure. The pre- 
sence of a heated chip below reduces the heat transfer rate 
at the middle chip. This happens because of the heating of 
air before it comes in contact with the central chip. The 
lower chip acts as a disturbance for rhe upv/ard flow of air, 
but since the natural convection velociuies are very small 
this does not distort the flow sufficiently to have any 
augmentation in heat transfer for the central chip. The 
effect of pre— heating continues to dominate for the range 
of Rayleigh numbers studied here, resulting in a lowering of 
its Nusselt number. Figure 4.17 gives the variation of 
Nusselt number for the two heaters. Figure 4,18 gives 
the rise in temperature for all the three heaters for this 
configuration. Maximum temperature rise occurs over the 
central chip, Tf we compare the temperature rise curves 
for the heater 1 in C-H-C and C-H-H configuration, it can 
be seen that rise in temperature is greater in case of 
C-H-H configuration. This is due to the presence of two 
heaters below the passive heater. 

These experiments show that on a circuit board, the 
components which are sensitive to temperature should be 
mounted as close to the base as possible. Powered compo- 
nents in series should be avoided even i£ they have low 
heat dissipation^ 
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The presence of a heated lower chip not only reduces 
the Wusselt number over the middle heater, but also flatens 
the Hu-Ra curve. Hence it is preferable to operate chip 2 
at a lower heat flux because for any increase in heat flux 
the corresponding increase in Nu will not be significant. 

The lo^/er chip does shov/ a gradual increase in Nu with Ra 
and hence can be operated even at higher heat flxixes. The 
change in the position of the neighbouring passive plate does 
not have any effect on Nusselt numbers reoorted here. 

All the Three Chips Heated 

Energy is supplied to all the three chips in this con- 
figuration. Experiments have been done for a wall spacing 
of 3H, 

Figure 4.i9 shows the variation of Nusselt number for 
all the three heaters for this configuration. The maximum 
Nusselt number occurs over the lowest chip as in the case of 
C-H-H configuration. This is again due to availability of 
cold ambient air for the lower chip. In this case also the 
middle dtiip has a lower value of Nusselt number and a compar. 
tively flat Nu-Ra curve. The upper chip is found to have tb 
same Nusselt numbers as the middle chip. This is because ai 
rising from the lower chip passes over the middle chip, 
accelerates resulting in the velocity increasing to a large 
extent. Hence there is an increase in the heat transfer 
coefficient over the upper chip. This increase in velocity 
compensates for the rise in temperature of air and the 
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Nusselt number for rhe upper hearer does not fall further. 

'■Je observe that a drop rn heat transfer rate is not 
always necessary for a ccxnponent if a few powered compo- 
nents are mounted belcn^r it. This gives an idea about the 
olacemenfc of IC chips dissipating different amounts of 
thermal energy on a circuit board. The most crucial 
component# highly sensitive to temperature rise# can be 
mounted on the lovjer side. In electronic systems there are 
always some components present \/hich dissipate a large 
amount of energy but can tolerate a reasonably high rise in 
temperature. Such components can be mounted at the top as 
the energy released by components below will not have any 
significant effect on them, further large temperature rise 
for such components v/ill not deteriorate the performance of 
the sysrem. 

Comparison with Existing Work 

Afrid and Zebib (1387) have carried out a numerical 
analysis for a two dimensional single and multiple uniformly 
heated devices mounted on an insulated wall. A two dimen- 
sional# conjugate# laminar flow model has been used and 
solid and fluid physical properties are assumed constant. 
They have reported the results for a heater with AR = 2. 

The Rayleigh number defined in their study is# 
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Ra = g ^ ar ) 
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Fig. 4, 21 : Configuration used by Afrid and Zebib 


If we define Rayleigh number as mentioned above the 

range of Ra for a single heater with AR = 2 in present 

6 6 

case varies from 0.25x10 to 0.52x10 • The corresponding 

variation of Ku is from 10.75 to 18.19. The value of Nu 

obtained by Afrid and Zebib varies from 26 to 42 for a 

6 6 

variation of Ra from 4x10 to 33x10 . In present case 
the value of Ra is limited by the temperature rise and 
power drawn by the heater. If we fit a curve# Hu = cRa^ and 
and calculate the constants c and ra for their data the value 
Nu for our Ra range varies from l5 to 16. In present e:5q)erime 
study we have obtained a larger increase in Nu with Ra. The 
reason might be that the slope of Nu-Ra curve decreases with 
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increase in Ra and the e:<ponents obtained on the basis of 
higher R?' '.ji-ll give c'^aller slooe. 

Afrid and zebib also shov; the effect of spacing between 
the components and the role of non~pov;ered and highly povjerec 
components. They have noticed that V7hen all the three 
components are heated the temperature rise between upp^r two 
components is smaller than the rise between the lower and 
the middle component. These results have also been obtained 
in the present experimental study. Experimental vjork by 
Orrega and Moffat (1985) also confirms thxs result. They 
report no change of temperature of the device after the 
fourth row of a set of 10 rows of aluminium components. 
Wximerical analysis of Afrxd and Zebib shows that the temper- 
ature rise varies almost linearly %rith the heat input for a 
single heater. This is exactly what has been found in the 
present work. 
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Forced Convection Experiments 
Single Chip 

Here we present results for a single heater with aspect 
ratio AR = 3, mounted on the wall of the flow apparatus. The 
range of Reynolds number covered is from 0 -• 1050. This 
range of Reynolds number has been studied because the veloci- 
ties encountered in the electronic equipment are small even in 
the presence of a fan. The range of Rayleigh number covered 
is 3000-65000. 

In the forced flow experiments different regimes such as 

forced convection# mixed convection and strong mixed convection 

regimes are encountered. These regimes can be distinguished 

2 

on the basis of the parameter Ra/Re , On the basis of the 
value of Ra/Re^ we classify^ 

2 

Ra/Re 0.1 as forced convection# 

^ , 2 

0.1 ^ Ra./Re ^ 1 as mixed convection# 

2 

and Ra/Re 1 as strong mixed convection. 

Heat transfer rate increases substantially in the presence 
of an imposed flow. Fig.4.22#4,23#4.24 and 4.25 show the vari- 
ation of Nusselt number with Reynolds number and different 
flow regimes at varying levels of input heat flux. 

At a low heating level# the Rayleigh ntiraber is small 
and buoyancy effects are insignificant . The heater goes into 
the forced convection regime# just after initiation of flow. 
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The slope of Nu-Re curve in this region is high and the 

increase in Nu with Re is large. As the input heat flux 

increases the Nu-Re curve becomes flat reducing the increase 

in Nu with Re. when heating levels are sufficiently high 

buoyancy effects become significant. Initial value of Ra at 
2 

q — 532 .2 w/m is 34517.4. This is the point bevond which 
buoyancy seems to augment the flow and hence heat transfer* 
Figure 4*26 shows the comparison Nu-Re variation at different 
heat fluxes. 

If we compare the variation of Nu at q = 68,97 and 
2 

q = 532,2 w/m , it can be seen that beyond Re = 200 the value 

2 

c£ Nu for q = 60,97 w/m is more than the value of Nu for 

2 2 
q = 532,2 w/m , This is because at q = 68.97 W/m , heater is 

mostly in the forced convection regime and the slope of the 

2 

curve is high. At q = 532,2 w/m the heater is in mixed 
convection having a flat Nu-Re curve. If the curves are com- 
pared in the same flow regime# the magnitude of Nu will be 
seen to be larger at higher heat inputs. This is seen in 
values below Re = 200 when both the heaters are nearly in 
same flow regime. 

At q higher than 532.2 W/m the buoyancy effects are 
important and we observe an increase in Nu with increasing 
heat flux. This can be seen fron the curves for q = 532,2, 
1006.6 and 1852.48 W/m^. 
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Tnis behaviour is different from what is known to occur 
in flush mounted heaters. The phenomena of heat transfer from 
a protruded heater is complicated. In such a flow field two 
regions of separation exist. The first region starts from 
the upstream wall and reattaches at the leading edge of 
heater. The second region starts from the trailing edge 
and reattaches on the downstream side of the wall* Heat 
transfer characteristics depend very much on the behaviour of 
these recirculation regions. These are affected by relative 
magnitudes of Rayleigh and Reynolds number. At low heating 
levels the buoyancy effects are small. The eddy size in the 
recirculation zone beyond the heater is small and weak. As 
the Reynolds number increases the size and strength of eddies 
on the upstream and down stream side of the heater increases. 
The large eddy gives better mixing with the colder fluid 
resulting in a higher rates of heat transfer. The flow in 
the upstream side separates from the wall and reattaches on 
the leading edge of heater. The leading edge is subjected to 
a colder fluid because of better heat transfer in the recir- 
culation zone before the heater. The thermal boundary layer 
thickness starts from zero at the leading edge. As the 
Reynolds number increases, these phenomenon become stronger 
and stronger resulting in a higher heat transfer rates. 

At a high heat input the initial Rayleigh number is 
large and buoyancy effects are important. The buoyancy force 
suppresses eddies in the upstream and downstream side. The 
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reattachment length on the downstream side reduces. The 
strength of vortex beyond the heater increases but because 
of a smaller size it results in poor mixing x^ith the cold 
fluid away from the heater .When reverse flow occurs at the 
heated supporting wall) energy is transferred to the main flow. 
This results in elevated temperatures in the recirculation zone 
and a reduced heat transfer coefficient. The eddy on the 
upstream side is also suppressed. Hot air gets trapped just 
below the heater. This heated film of air is stably stratifi- 
ed and reduces heat transfer. The approaching cold fluid is 
also affected by the trapped heated air and its temperature 
increases. The leading edge of the heater is now subjected to 
heated air and the thermal boundary layer over the heater 
does not start from zero as initial thickness. This results 
in reduced heat transfer over the vertical front surface of 
heater. The buoyancy force does tend to increase the velocity 
over the heater but the net augmentation in the heat transfer is 
less. A larger thermal boundary-layer thickness results in 
small temperature gradients normal to the mean flow direction 
and a weakening of the buoyancy force hence the heat transfer 
is less than that obtained in forced convection regime. This 
is the main reason due to which a higher Nusselt number is 
obtained in the forced and weakly mixed convection regime. 

As the input heat flux is increased beyond a certain 
point buoyancy effects begin to dominate and improvement in 
Nusselt number is seen. This explains the increase in Nu when 
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, 2 

the heater xs xn a hxgher Ra/Re range for q = 1006,6 and 
q = 1852,4 w/m , 

Sxmilar behaviour has been reported by Habchi and Acharya 

(1985) in their numerical work on laminar mixed convection in 

a partially blocked# vertical channel. They have investigated 

the flow and temperature field numerically at different Rayleigh 

2 

numbers for different values of Ra/Re , They report that at 
a given Ra# the recirculation region above the blockage becomes 
stronger and smaller as Ra/Re^ increases. At a high Rayleigh 
nxamber# Nu along the upper surface of the blockage increases 
as Ra/Re^ increases due to an increase in the recirculating 

3 

eddy strength. At a lower Rayleigh number (Ra = 10 ) they 

2 

observe a decrease in Nu along the upper surface as Ra/Re 
is increased because of very small and weak recirculating eddy 
beyond the blockage. They have also reported the values of 
average Nusselt number. It is found to increase with increas- 

y 2 

ing Ra/Re at high Rayleigh numbers. At low Ra and Re# the 
recirculation zone is very small and its effect is not seen 
on average Nusselt number. 

Shou-Shing Hsieh# Huei-Jan Shih and Ying-Jong Hong (1988) 
have presented the results for laminar forced convection from 
surface mounted ribs in a horizontal channel. They numerically 
study the low speed forced convection heat transfer near two 
dimensional transverse ribs in which the walls are held at a 
uniform temperature. 

They define the average Nusselt and the Reynolds number 
on the basis of heater depth H. They have also considered 
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for the cases when mixed convection effects are negligible 
(at q = 40.64/ 68.97/ 102.05 and 126.96 w/m^) the Nusselt 
number is found to vary within the range of 6,89 at Re^^ = 84 
to nearly 14 at Re^^ = 333. The range of Nusselt number cal- 
culated experimentally in our case is found to agree well 
with the results presented by Hsieh et al, despite the present 
study being applicable for a vertical configuration. The thermal 
boundary layer thickness over the heater in case of forced conve- 
vection is small than in natural convection. j^nj increase 

in the distance between the he-ater and the passive plate beyond 
2H does not affect the heat transfer. 

Since Nu increases with increasing Re and the increase in 
Nu is greater in forced convection regime/ it is preferable to 
operate with a large Re to have good cooling performance rather 
than rely on augmentation due to mixed convection. 

Effect of a Protrusion on the Opposite Wall 

Metallic strip 


Pig.4.27 1 

We have studied the effect of a baffle mounted on the 
opposite wall in case of forced convection experiments with 
single heater. We present data for the three values of heat 
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fluxes (q = 104,43, q = 222.5, q = 783.77 W/m^ ) . Vlith a 
sharp protrusion in the way of flow the resistance to flow 
increases. Hence the maximum velocity in the test cell is 
reduced limiting the largest Reynolds number to 780, 

Fro?! Figure 4,28 it can be seen that the buoyancy aug~ 
ments the heat transfer more clearly with a baffle. No 
improvement in the Nusselt number is observed at q = 104,43 w/m 
if the results are compared with those for a heater without 
protrusion on the opposite wall. This can be seen from the 
Figure 4,29 when the two ^^a.ie.Ys are at nearly equal heat 
input. Figures 4,30 and 4,31 show the comparison between 
the two cases at higher heat flux, (Note that there is some di 
fference in the energy input between a flush and a protru- 
ded heater) . 

At higher heat fluxes the Nusselt number with and without 
a baffle is found to be nearly equal. If we look at the 
curves the value of q for heater with protrusion is 222,5 
w/m and for heater without protrusion on opposite wall 
is 532,2 w/ra , Because of reduced buoyancy effects the 
expected Nu-Re curve at q = 222,5 W/m in case of heater 
without protrusion is likely to fall on the curve for the 
heater with a protrusion on the opposite wall. By a 

2 

similar reasoning the expected curve for q = 783-»77 W/m for 
heater without protrusion will fall almost on the curve 
for heater with protrusion. 




for the cases with and without protrusion on opposite wall 
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The absence of augmentation in heat transfer in the 
presence of a protrusion on the opposite wall can be explained 


as follows. 


Recirculation zone 
'on do\^m stream side 


Recirculation zone 
on upstream side 





Separation behind strip 


metallic strip 


Fig, 4, 32: Flow pattern with protrusion on opposite wall 
Fluid rising in the channel is diverted towards the heater in 
the presence of the protrusion on the opposite wall. This 
diversion of fluid is likely to have suppressed the vortex 
in recirculation zone on the upstream and downstream side of 
heater. Though the velocity near the heater increases because 
of reduced area, the lower heat transfer from the upper and 
lower sides of heater may be the reason for no augmentation 
in Nusselt number. 


Since the fluid flowing is closely packed near the 
vertical surface buoyancy seems to have a stronger 

effect in this case shovjing a clearer rise in Nu with 
increased heat flux. The recirculation regions being 
already suppressed in case of heater with protrusion an 
increase in the slope of Nu-Re curve is not seen at low 
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heating levels. 

The placement and the size of baffle will have an effect 
on heat transfeir frcsn the heater. The presence of protrusion 
increases the resistance to flow substantially. Hence it may 
not be really economical unless there is marked improvement 
in Nusselt number. This aspect must be investigated in 
detail in future studies. 

Multichip Forced Convection Experiments 

In these experiments forced convection studies over a 
multichip heater have been carried out. The configurations 
considered are: 

Central chip heated 
Upper and lower chip heated 
All the three chips heated 
Upper two chips heated. 

In these experiments, due to increased resistance to the 
flow and the reduced size of heater, the range of Reynolds 
number is reduced to 550. The variation of Rayleigh number 
is up to 15000. Maximum temperature encountered is nearly 
100*^C, Foxir different heat flux levels have been considered. 
The aspect ratio for these heaters is 2, 

Central Chip Heated 

Figure 4.33 shows the variation of the Nusselt number 
with Reynolds number at different heat fluxes for the central 
heater. 



At q = 181*3 w/m the curve is initially flat and then 
rises as the Reynolds number is increased. The other curves 
are in the mixed convection regime having a flatter response 
With Re* The range of Re is from 0 to 550 only. At 
these low Reynolds numbers buoyancy will have a stronger 
effect. It is also due to differences in the flow pattern 
in comparison to a single chip heater. In the valley between 
the lower and the middle chip a recirculation zone exists 
because of flc«,«; separation over the lower chip. This 



Fig.4.34 i Flow patberit—for nultichip heater 
recirculation zone will lead to a better mixing of fluid in 
this region and an irt^jroved heat transfer* This is to be 
contrasted with a single chip heater where a hot stagnant 
film is formed downstream* Similarly the vortex in the 
region between the upper and middle chip is also not suppressed 
by increase in buoyant force as its strength primarily depends 
on the strength of flow. At higher Reynolds number the mixed 
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convection effects become insignificant and all the curves 
tend to merge. 

Upper and Lower Chip Heated 

The problem has been considered for 4 different heating 
levels. It is seen that the Nusselt number for the heater 1 
IS more than the Nusselt number for heater 3 and the differ- 
ence reduces at higher heating levels , This can be seen from 
Figs. 4.35 and 4,36.The lower and middle chip disturb the flow 
and an increase in heat transfer over the upper heater is 
obtained. At higher heating levels the temperature of fluid 
goes up when it passes over the lower heater and hence the 
difference between the Nusselt number for the two heaters is 
reduced. 

Fiq.4.37 and 4.38 present the variation of Nu with Re at 
various heating levels for heater 1 and 3. Since the Re range 
IS small, buoyancy effects are clearly seen. For heater 1, 
the augmentation due to buoyancy decreases as the heat flux 
increases. At a high input heat flux the temperature of 
fluid passing over the heater 1 is sufficiently high and 
uniform because of heat gain from the lower heater. This 
reduces the temperature gradients in the direction perpendi- 
cular to the heater. Subsequently there is reduction in 
augmentation due to buoyancy. 
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All the Chips are Seated 

we have studied this configuration for three heating 
levels • Pigs. 4. 39 to 4.41 give the comparison of 

Nusselt number for the three heaters at different heat flux. 

In all the cases it is seen that the heater 2 has the lowest 
Wusselt number. Heater 3 is found to have the highest Nusselt 
number as it is subjected to the cold ambient fluid. The 
lower heater increases the disturbance in the flow over the 
middle heater but the temperature of fluid also increases. 

Since the value of Re is low the effect of rise in temper- 
ature of fluid is predominant and the Nusselt number over the 
middle heater drops a little. As the fluid passes over the 
central heater the disturbance to the flow increases signi- 
ficantly resulting in a better heat transfer over the heater 1. 

Nusselt number of heater 1 is nearly equal to that of heater 3 

/ 2 

at q = 192,7 and 634,5 W/m , At a further increase in heat fit 
the temperature effects are important for heater 1 resulting 
in Its Nusselt number being equal to that of heater 2. 

Figures 4.42 to 4,44 show the Nu-Re variation indivi- 
dually for the three heaters at different heating levels. Aug- 
mentation due to buoyancy can be seen for all the three heaters, 
As the Reynolds number increases, the forced convection 
becomes significant and the curves seem to merge. For the 
upper heater, as in case of H-c-H configuration, the effect 
of buoyancy seems to reduce because of the reasons given 


earlier 
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Fig* 4*43 





McEntire and webb (1388) have studied experimentally 
forced convective heat transfer from protruding and flush- 
mounted two dimensional discrete heat sources. They have 
measured local heat transfer coefficients vjithin the range 
of 1000 ^ ^^e ^ 10000, For the protruded heaters they find 
enhanced heat transfer for the second and subsequent heaters. 
This is explained by the flow separation and associated venna 
contracta effects, with eventual downstream transition to 
turbulent flov; induced by the presence of protruding heater. 

Similar disturbing effects have been seen in the present 
study. Since in our case the range of Reynolds nxiraber is 
small the temperature effects are also significant and 
we do not see any transition to t^urbulence. This is the 
reason for slight lowering of heat transfer coefficient over 
the second heater (heater 2) . The heat transfer coefficient 
for the third heater (heater 1) is nearly equal to that for 
the first heater (heater 3) which shows that at higher Rey- 
nolds numbers the disturbance initiated by the first heater 
will increase and subsequently an enhancement in heat 
transfer coefficient for downstream heaters will be obtained. 

Upper Two Chios Heated 

From the Figure 4.45 it is clear that the Nusselt number 
for the heater 1 is more than that for heater 2, In this 
case the heating level is low and temperature effects are 
small. The flow passes over the lower two heaters and the 
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disturbance is high enough to enhance the heat transfer overt:e 
upper heater# Since we do nor have data for a higher heat 
flu:<, we can not present the exact behaviour in this config- 
uration but the difference between the Husselt numbers for 
two heaters is expected to reduce because of increasing 
temperature effects at higher heating levels. 

Let us consider the effect of a cold protrusion mounted 
before a heater. For this study we look at heater 2 in C-H-C 
combination. This heater is subjected to an unheated protru- 
sion and cold ambient flow. The Nusselt number variation 
for this heater can be compared with that for the heater 3 
in H-H-ri and H-c-H configuration. The heater 3 in these con- 
figurations is exposed to cold ambient flow directly. The 
value of Nu 2 for C-H-c configuration is in general higher 
than Nu 2 for H-H-H and H-C-H configurations if compared for 
nearly same heat flux and Reynolds number. 

This shows that the presence of a cold protrusion in 
the path of flow before the heater increases the distur- 
bance and hence heat transfer. This suggests that the placing of 

1 

an unpowered component before a powered component will lead to 
a better thermal performance of the electronic system. 

On comparing the value of Nu^ for H-C-H and H-H-H 
configuration, it can be seen that the placing of cold com- , 

ponent above a powered component does not have any significant 
effect on the heat transfer. 
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If we compare the Nu^ in H-c-H, H-H-H and H-h-c config- 
urations^ it is seen that maximum Nusselt number is obtained 
in case of H-C-H configuration in general. This gives a 
design iaea about the placement and spacing of components in 
an electronic system. In these configurations the heater is 
subjected to disturbance created by two protrusions on the upstream 
side. Nu^ in H-C-H is more than Nu^ in H-H-H because a cold 
chip is placed between the two heated chips reducing the 
heating effects of the lower chip. Wu^ in H-H-C is more 
than ^^^2, H— H— H as in the former 

the heater 1 is preceeded by one heated protrusion. In H-H-H 
heater 1 is preceeded by two heated protrusions. Heat transfer 
will increase if an unpov/ered component is mounted between two 
pox'/ered components. Alternatively the distance between the two 
IS increased which will result in a better mixing with an 
entrainment of cold ambient fluid. 

Let us compare H-H-H configuration in forced and natural 
convection. In both cases the maximum Nusselt nximber occurs 
over heater 3, In natural convection the difference between 
Nuj and Nu 2 is large and then the Nusselt number becomes nearly 
constant for further heaters located downstream. Owing to 
disturbance in the flow the difference between the Nu^ and 
Nu 2 is small in forced convection. The Nusselt number incre- 
ases further for the heater 1, It is expected to rise further 
for the heaters located further downstream. This signifies 
the complicated flow phencanena and role of disturbance in 
the flow initiated by the heaters mounted on the upstream side 
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in forced flo^TS. In both cases location, spacing and 
size of components are quite important. In forced conve- 
ction flow is complex over and between the heaters. Besi- 
des mixed convection effects, recirculation and separated 
regions exist and analysis and interpretation of data are 
more involved. 





Qualitative Analysis by Interferometry 

In this section we present a qualitative description of 
the thermal field over a single chip and multi-chip heater. 

When the interferometer is initially in the infinite fringe 
setting, each fringe is an isotherm. These isotherms give 
the idea about the thermal boundary layer. To a limited 
extent flow pattern, recirculation and flow separation can also 
be predicted. In the present study each isotherm corres- 
ponds to a change in temperature of 

For a single chip heater all the photographs have been 
taken from the monitor of image processing system with a 
highly sophisticated camera (Canon T90) with a 400 ASA film. 

Some of the pictures are pseudo-coloured to show a clear 
contrast between fringes. Another advantage of pseudo- 
colouring is that different temperature zones can have different 
colours and thus one can have a clear observation of the 
temperature field. In present study we have not made use of 
this facility and the colouring has been done only to make 
the picture sharp* 

For the multichip heater, the photographs have been taken 
by the field camera and the film is developed in the labora- 
tory itself. This causes some of the pictures to be not very 
sharp. Moreover thermocouple wires in the path of the light 
beam also distort the image by producing their own dispersive 
optical effects* 


d 
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Single Chip Heater 

Both the transient and the steady temperature field 
inside the enclosure have been visualized using interfero- 
metry. Figures 4,46 & 4,47 show the isotherm over the 
heater during the transient and steady state respectively. 

The first film has been taken 20 minutes after the cold 
initial state. The number of fringes is fewer and the 
fringe spacing is large. This shows low temperature gra- 
dients and hence a low Nusselt number during the transient. 
The number of fringes at steady state is fairly large. 

These are closely packed, indicating a large temperature 
gradient and a high Nusselt number. Thermal boundary layer 
can be seen in the Figure 2. The thickness of the thermal 
boundary layer is small at the leading edge of chip and 
increases over the heater. This shows that the highest 
value of the Nusselt number occurs at the leading edge. 

The air outside the boundary layer is cold as seen by the 
absence of fringes in that region. 

Figures 4,48 and 4.49 show the trailing side of heater 
for a low and a high input heat flux respectively. The 
thickness of fringes increase near the trailing edge due to 
a reduction in heat transfer in that region. On the down- 
stream side of heater, the fringe spacing is increased which 
is a region of low velocity and low heat transfer. The 
deviation of the outer fringe indicates the existence of a 



recirculation zone. Further downstream the deviation of 
the outer fringes towards the passive wall is also observed. 
The isorherms remain attached to the passive wall beyond 
the trailing edge. Hence w'e expect the temperature in 
this region to be higher. If cooling is provided here 
then the thermal performance of the circuit board can cer- 
tainly improve. 

Figure 4.50 shows the isotherms in the upper region 
of the enclostire when the top of the enclosure is open, 
rfere a v/all plume like behaviour is seen and the thermal 
energy is lost to the ambient from the top instead of 
accumulating on the side wall as in the case of a closed 
enclos ure. 

Figure 4.51 shows the initial cold state with the 
interferometer set in the wedge fringe mode. Subsequent 
deviations of the isotherms at steady state is shown in 
Figure 4.52 . Here again a boundary-layer can be seen 
whose thickness is equal to that in the infinite fringe 
setting for a given input heat flux. 




Isotherms in steady state 






Fig, 4,48 : Isotherms downstream the heater at low 
heat flux 



T- 14 „ /I /in . f^nwn<=!t"np«Tp the heater at high heat flux 


^i9«4,50 ; Isorherms in the upper region of an open 
enclosure 






^ig* 4.51: Initial setting of interferometer in wedge fringe 
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Fig* 4*52 i Wedge fringes at steady state 


Mulichip 


The proble.n can oe divided into 3 configurations. 

(a) Central chii* hratcd 

(b) Letter t\;o chips <-re heated 

(,c) All the three chips are heated. 

Central Chip Heated 

/ 2 

The heat flux applied is q = 7J9 W/m and temper- 
atures are T^^ = 53 °C and T^ * 27 °C, The transient state 
can be seen in figures 4.53 and 4.54 . In the initial 

stages the isotherms tend to tilt away from the heated 
strip because of weak flovj. As the temperature of chip 
increases rhe buoyancy force increases and the natural 
convective flow in the vicinity of chip becomes more para- 
llel. 

Figure 4.55 shoves the middle chip at the steady state. 
The growth of thermal boundary layer in the vertical dire- 
ction can be seen here. 

Figure 4.56 shows the lower chip and the spacing between 
middle and lower chip. Since -the lower chip is unheated only 
one fringe ic formed on its surface. Most of the fringes can 
be seen originating from the wall which shows the heating of 
wall close to the chip. The absence of fringes in the lower 
portion of enclosure indicates a cooler region os .ejgjected. 




Figure 4.57 shows the isotherms in the region between 
upper and middle chip. The isotherms very close to the 
heated chip merge in the region close to the trailing edge 
of the chip. This shows that the temperature of chip is 
much higher than temperature of the supporting wall which indi- 
cates that conductive heat transfer to the bakelite sheet is 
low. The upper strip is also unheated except by the pltome 
ccxning from below* Hence it does not have many isotherms. 


_ j J'iC r \c OV'- 


^ig.4.54 

r central chip in tr^n- i^-nt - tar 


Fig. 4. 53 




Fig, 4. 55 : Isotherm^ ct £-~eady 
state 


Fig,4,5C : Spacing benvjeen middle 
and lox^er C^xp 


Pig. 4.57 


Isotherms in spacing betv^en upper and middle chip 


Pne prccaniic o le-’ is-, -her. is bovever does show the tem- 

veroture oZ u > ’ c' :: ir hi her than the ambient temper- 
ature. Sin'-c UMier chip does not dissipate any heat flux 

the isotherms ‘-C- from the chip. No significant 

thermal boun-.ar-- l^^yar -rrowth is seen. The number of isotherms 
in the lov’cr r.paon of enclosure being less than that in the upper 
region cb'arly in.'i'-ates the difference between the heat 
transfer mechanisms to loi;er and upper heater and the dire- 
ction Oi. mean buoyant flow. The lower chip gets heated because 
of diffusion only and it is reasonable to expect the wall 
con,<uction to be insxcniiicant . 


Lo\;er Tx;o Chios are Seated 

in this configuration the middle and the lower chip are 
dissipatincj heat. The chips are subjected to only moderate 
haar flue levels so that the fringes formed on them can be 
seen clearly. The figures 4.58 to 4,61 show the 
thermal fielu close to the lower chips in transient state, 
jj'igure ^,5.' corresponds to initial infinite fringe setting of 
interferometer. Figure 4.59 shows the beginning of the for- 
mation of isotherms. This film has been taken after 50 minutes 


from initial cold 

at steady state. 

central chip is <32 

q3 = 208.4 • 

= 23,8°C, '3?2 = 


state. Figure 4.61 shows the fringes 

In this configuration the heat flxix on 
2 

= 202.6 w/m and on lower chip is 
34.72°c, T3 = 31.96°c and = 


18.3°C, 




Fin.4.5S:Initial infinite Fig. 4.59: Tranrt~^_^ ^ 

The isotherm closest to the heater follows the wall 
along the heater surface. It follows the wall in 

the valley also between the two chips. Since in this case 
the heat fltix is low, the flow is weak and separates from 
the heater surface. Higher order fringes show early separa- 
tion by tilting away from the wall. 



-I'". '.CP: rrcmsicnt state Fig. 4. 61 : :t:cady state 

It can be seen frcm the Figure 4.61 that rate of heat trans- 
fer and subsequently the Nusselt number is highest at the 
leading edge of the lower chip. The central chip does not 
come in direct contact with cold air. This is seen from 
the isotherm that begins from the lower chip and envelopes 
the middle chip. This explains the lowering of Nusselt 
number over the central chip. 

The Figure 4.62 shows the isotherms in the region 
between middle heated and upper unheated chip. Since the 
upper chip is unheated/ separation occurs right at the 
trailing edge of the central heater. It can be exolained 
as follows. Consider the isotherm which follows the wall 
and the valley between the two heated chips. It begins to 
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361'urate at the end of central heater. 

The ij'xgure 4,t53 shows isorherms in the upper region 
of enclosure, ag the upper chip is unheated the isotherms 
snou a deviation towards the neighbouring wall and accumu- 
lation of hot air in the upper region. It can be seen that 
the uiper cold chip acts as a disturbance to the thermal 
plume rising fron the central chip resulting in the iso- 
therms get deviated towards the opposite wall. This energy 
can be absorbed efficiently ij the cooling and ventilation are 
provided on the opposite v/all* 

All The Three Chios are Heated 

In this case the chips were subjected to a higher heat 
flux. As can be seen from the figures that the fringe density 
IS very high close to the chips. As in the earlier config- 
urations the maximum Nusselt number occurs at the leading 
edge of the lower chip since the thermal boundary layer 
thickness is a minimum at this point. The heat flux and 
temperature levels are- 

2 

= 1334.6 w/m^; q.2 = 1319.0 w/m 
= 1300.5 w/m^ 

T^ = 87.02°c ; Tg = 86.8°C ; T3 = 64.16°C 

Too = 

Figure 4.64 shows the isotherms over the lower and central 
chip. The region between the two chips also becomes hot. 



order isotherms close to the heater follow the 
contour Of the heater and the valley between heaters. 

ws the region between the central 
and upper heater. It can be seen from this figure that 
the temperature of air inside the enclosure increases as 
3V upwards. Figure 4,66 shovjs the isotherms in the 
Upper portion of the enclosure. The separation of isotherms 
rrom the wall is seen here. The deviation of isotherms and 
congestion in the upper region indicates hot air accumulation 
and some amount of air-circulation. 


1 he lower region of enclosure is shovm in Figure 4.67 • 

It can be observed that the temperature of the enclosure in 
this region has not risen much above the ambient temperature. 
It is only beyond the middle of the central heater that the 

mean temperature of air inside the enclosure rises subs- 
tantially. 


As we begin to move up in the enclosure, it can be seen 
that the higher order fringes follow the wall upto a point 
3nd then separate and deviates towards the neighbouring 
wall. Their thickness also increases along the path. The 
fringe pattern near the opposite wall indicates some amount 
of recirculation in the enclosure as the hot air on the active 
wall rises up and pushes the air trapped in the upper portion. 




The thermal boundary layer formed over the central and 
heater ere seen to be nearly equal in thickness. 

Hence the heat transfer rates and Nusselt number over the 
two heaters are nearly equal. This is in conformance 
with the result obtained in the energy balance experiments. 
There is a variation in Nusselt number over the heater 
length and this can be seen from the variation in the fringe 
density over the heater length. 


Fig. 4. 64 ; Isotherms between lower Fig. 4.65 
two chips 


Isotherms over upper 
heater 



uS have calculated the value of Musselt number by 
xf Lcr .'eronetr’, for a few cases to cross check the results 
ccoo-rnej. using energy balance. Here present results 
.or t o cases. The method used to calculate the Nusselt 
nai.ioer rs based on curve fitting referred earlier in 

w* a ^ o X. a . 


Case 


fhe neater asoect ratio is 3, 


Inr-ut hear _"lxr;c 


564.72 w/m 


Temperature of chip = 48.7 C 

o 

Ambient temperature 18.1 C 


?he Mucsolt number through conventional method = 18. 1. 

Ihe ITuesf'lt number through interferometry = 20.71. 


Case 2: 


T = 53°C 

c 

T = 27°C 


The heciter aspBcc rstio is 2© 

Input heac flux q = 799 w/m 

Cnip temperature - 53 C 

Ambient temperature = 27 C 

The Nusselt number through conventional method = 14.3. 

The Husselt number through interferometry = 15-1. 

The ilusselt numbers obtained by the two methods are close 
but the value obtained using Interferometry is higher. 




n j.r 


3r w r'ction length being small the estimated 
:oL ns ' li.',;iblc at: the chip surtace. This is 
* tc "f'Crature gradients near the heated xrall, 

Si.'*'* jr'client: reuuces away from the chip surface, 
'i-nrr coo'-en to determine the teifiperature profile 
-r > _r''i Lhf chip surfaces the temperature gradients 

' the error due ;o refraction at those points 
* rt*j.l t-n’ c^n be neglected. The accuracy of this method 
sHiiS a on the number of points chosen and their distance 
.'rom the hr.atGf surface. More the number of points better 
’.;ill oe tr/' fitted curve. The points taken near the heated 
rurface ’•”-11 give a better approximation to the actual tera- 
i ereturo cofile in the vicinity of the heater. Measurement of 
fringe s I'-cing manually through the negative of a photo- 
gre^phic film also involves some amount of error and approxi- 
mation. 

Considering all these factors/ the results obtained by 
the tv 70 methods can be considered to be close v/ithin experi- 
mental limits. 


Calculation of the Nusselt Humber (Nuiae gicall. 

In the Table 4-53 presenr the results for a 


flush heater mounted on the wall of a vertical channel as 
calculated numerically. The basic scheme is described in 
Muralidhar and Kulacki (1987). -Were the results have been 
obtained with some modifications in the boundary conditions 
The data presented gives an idea about the Nusselt numbers 
with a flush heater in mixed and forced convection regime. 


142 


.ret. - vjn\ w-,on hac nor been considered in numerical calGu~ 
^t.rion. - o- configuration and detailed numerical anal3''sis 
rc-ar A_ 4 . 

-- can toe observed that the ITusselt numbers \rxth a 
-lash n~c.cer are lox«;er than that with a protruded heater 
(uG Ccilculc-ced coiperimentally in this study)* This result 
is »sctc.d ioecause of better heat transfer rates obtained 
\;itn - crorruced heater. It is clear that as the Rayleigh 
number increases the Nusselt number increases for the s ame 
Reynolds number. Plow patterns are simple and the buoyancy 
has a clear augmentating effect on the Nusselr number for 
a flush heater. Some amount of recirculation on the opposite 
wall has been observed at high Rayleigh numbers. Experiments 
ith protruded heater show a different behaviour because of 
conplex mecnanisms of flow and heat transfer. 

Huraerical results have been obtained for flat and 
parabolic velocity profile in the channel. The Musselt 
number obtained with a flat velocity profile is liutle 
higher than with a parabolic velocity profile. But the 
velocity profile does nor affect the Nusselt number too 
significantly. 
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Table 4,1 : Single Chip Natural Convection 


Open Enclosure 


q 

T„ 

T 

Ra 

Nu 

c 

oo 


60.99 

35,4 

30.1 

8333 

8.732 

208.27 

41,38 

30,2 

14687 

16.2 


411.55 

51,95 

30.2 

24992 

17.2 


755.04 

53,6 

30.0 

34911 

20.45 


1178,46 

79.2 

30.1 

44663 

21.6 


1396,5 

90.89 

30.3 

49993 

20.75 

,D 


1752.9 

104.1 

30,5 

54547 

21.4 



D = 3H , AR = 3 


Table 4.2 ; Single Chip Natural Convection 

Closed Enclosure 


q 



Ra Nu 


63.49 

35.17 

212.9 

43.15 

432.56 

52,88 

775.08 

66.65 

1176.71 

83,60 

1388.35 

91.84 


31.0 

5458.9 

14.2 



29.6 

16904.2 

14,61 


33.0 

20313 .5 

19.66 

1 

30.3 

36728.4 

19.1 

J 

29.3 

47828.7 

19.0 


29,1 

51702.4 

19.0 



D 


2H, AR = 3 
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Table 4.3 


Single Chip Natural Convection 
Closed Enclosure 


q 

T 

T^ 

Oo 

Ra 

Nu 

52 ♦ 6 

35. 3 

30.1 

8208 

9.4 

■ 

203.4 

46.5 

30.3 

19529 

11.8 

1 

451.9 

55.09 

33.3 

23926 

18.75 

il 

774.3 

67.40 

30.2 

37367 

18.6 

■ 

1211.0 

84.90 

30,6 

46982 

19.5 


1390.6 

93.75 

30,2 

51258 

18,9 

- 

D = 3H, 

AR = 3. 





Table 4.4 : Single Chip Natural Convection 

With a protrusion mounted in the path of fluid. Location 
of protrusion is shown in Figure ( Ao, ) 


Position of 
disturbing 
strip case 

q 

^c 

T 

OO 

Ra 

Nu 

a 

440.45 

56.76 

33.83 

24751 

17.33 

b 

435.45 

58.49 

34.26 

256S7.1 

16.11 

c 

439.3 

58.42 

31.52 

28982,9 

14.74 

d 

439.3 

58.31 

31.31 

29153.11 

14.74 

e 

437.0 

57.92 

31.8 

28220.3 

15.15 

f 

434.32 

58.30 

31.14 

29357,1 

14.47 

g 

435.44 

58.12 

31.0 

29383.4 

14.53 

h 

437.07 

57.33 

30.5 

29354.2 

14.75 

i 

1169.3 

86,32 

31.9 

46211.0 

18.59 

j 

1173.8 

86.95 

32.5 

45848.6 

18.65 

k 

1169,3 

87.01 

32.1 

46313 .6 

18.48 


AR 
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Table : Single Chip Natural Convection 


White Enclosure 


q 

T 

c 

T 

oo 

Ra 

Nu 

337o9 

50.27 

29.95 

23799.8 

17.5 

716.6 

55.55 

30.9 

35223,5 

18,55 

1153.4 

82.43 

31.81 

44407.2 

19.91 

1487.5 

95.21 

33,3 

48536,8 

20.67 J 

1432.4 

91.41 

29.51 

51075.9 

2 0 # 82 £> 

1497.8 

94.82 

31.3 

50493.3 

20.36 

1333.45 

104.66 

29.85 

55388.3 

20.965 

2213.77 

115.98 

30.26 

57945.8 

21.45 


D = 3H, AR = 3 


Table 4.6 ; Single Chip Natural Convection 

Black Enclosure 


179.3 38.73 

202.8 38.22 

341.15 46.43 

719.47 67.78 

1167.31 80.39 

1509.7 95.62 

1928.99 106.19 


T 

Ra 

Nu 

oo 



29.5 

11993.6 

17.94 

28.1 

13321.7 

18.3 

29.25 

20851.7 

13.35 

33.87 

33261,5 

18.95 

34.11 

40776.2 

22.104 

34.47 

47485.8 

21.20 

33.93 

51733.6 

22.62 









Table 4.7 i Single Chip Natural Convection 
Radiative Heat Loss 


Chip 

temp. 

0-. 

C 

Power 

input 

watts 

Radiation 

losses 

watts 

% Radiation 
losses 

38.73 

1,434 

0.1433 

9.99 

45.43 

2.72 

0.2797 

10.2 

57.7 

5.745 

0.6436 

11*83 

30.38 

9.321 

0.9748 

10 . 4'5 

95,2 

12.05 

1.533 

12.7 


D = 3H, AR = 3. 


Table 4.8 ; Single Chip Natural Convection 
Setup without Slit 

(with Conduction to supporting plate) 


q 

^c 

T 

oo 

Ra 

Nu 

Su 

AR 

42.9 

33.04 

30.37 

3599 

15.101 

5.03 

38.3 

35.43 

30,5 

6533.15 

15.75 

5.55 

185.6 

38.42 

28.65 

12821.6 

17.67 

5.89 

352.15 

46.67 

28.78 

22025.6 

17.99 

5.99 

685.8 

60.87 

29.1 

35049.0 

19.34 

6.44 

840.2 

70.27 

29.45 

41846.0 

18.21 

6.07 

1172.9 

79.9 

29.6 

48293.0 

20.39 

6.79 

1463.4 

89.4 

30.0 

53013.0 

21.285 

7.09 


Ra 


AR 


3 


133.2 

241.9 

474.8 

815.7 

1298.1 

1549.8 

1788.6 

1963 . 4 ^ 


D 


3H, AR = 3 



Table 4.9 : Setup with Slit 
(:io Conduction to the supporting plate) 


q 

c 

oo 

Ra 

Nu 

4 3.36 

33.81 

30.37 

3964.6 

14.81 

37.341 

36.72 

31.1 

7333.9 

16.28 

130.85 

37.57 

30.5 

9215.2 

17.25 

374.39 

49.26 

28.78 

24801.0 

16.67 

705.03 

67.21 

29.25 

39971.0 

16.54 

983,2 

77.53 

29.5 

49940.0 

17.96 

1227.45 

87.96 

29.6 

52948.0 

18.25 

1503.88 

97.00 

30.02 

57070.0 

19.23 


D = 3H, AR = 3. 


Table 4.10 : Single Chip Natural Convection 
Heater with Reduced Aspect Ratio 


^ "^chip 

141,6 35,99 

329.49 43.29 

791.3 58.23 

1453.0 77.29 

2135.4 98.75 

3004.3 118.93 

3859.6 133.8 


T 

OO 

Ra 

29.7 

2106.21 

28.9 

3888.7 

23.7 

7006,6 

28.7 

9776.3 

29.4 

11654.4 

29.9 

12729.7 

29.0 

13441.1 


Nu 

AT 

Hu 

AR 

Ra/AR 

10.75 

7.29 

5.37 

263 

12,74 

14.39 

6.37 

486 

14 .59 

30.93 

7.2? 

87 5 

15.890 

48.59 

7.942 

1222 

15.892 

69.35 

7.945 

1456 

16.96 

89.03 

8.48 

1591 

18.19 

104.8 

9.09 

1680 



D = 3H, AR = 2 


i-’aole 4,1 1 ; Natural Convection 3-D Heater 
oantcd on ahorisontal Plate .All sides exposed to ambient 


q 


oo 

Ra 

Nu 

m: 

2355 

53.55 

25.2 

4898 

32,2 

33,5 

2765 

52,3 

25.2 

5363 

33,2 

37.7 

3214 

67.5 

26,6 

5574 

35.32 

40,94 

34S4 

72.83 

23,4 

5585,9 

35.65 

43.49 

5323 

33.21 

23,79 

72 08,0 

35.87 

64.42 

5366 

33.52 

29,31 

7495 

36.5 

70.21 

5573 

107.14 

29,08 

7925,5 

35.84 

78.06 


Ilountcd 

Table 4.12: 

on a Vertical 

bJ 

Natxiral Convection 3-D Heater 

Plate, All side exposed to ambient 

q 

c 

T„ 

oo 

Ra 

Nu 

Zffi 

2195 

54.09 

26,2 

4202 

36.1 

27.9 

2812 

53.50 

25.8 

4893 

38.0 

33.7 

3 42 6 

66.89 

26,8 

5476 

38.47 

40.1 

3957 

76,25 

26,9 

6299 

35,6 

49.4 

4101 

78.24 

27.5 

6359 

35.8 

50.7 

5291 

84.61 

28.4 

6693 

41.2 

56.2 

6069 

91.66 

26.0 

7566 

40.34 

65,66 


rable .'.13 ; Natural Convection 3-D Heater 
Vertical Configuration Enclosure 
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Spacing 

D 

q 

c 

T 

oo 

Ra 

Nu 

LT 

°C 


1053 

45.6 

33.3 

1847 

39.46 

12.32 


1553 

51.7 

33.1 

2703 

40,59 

18.56 

5H 

2280 

59.9 

33.6 

3596 

39ol 

26.29 

2533 

52.21 

33.6 

3844 

40.65 

28.5 


3539.7 

67.5 

33.1 

4472 

45.9 

34.42 


52 57 

94.6 

33.2 

6720 

42.4__ 

61.4 


1354 

49.2 

32.7 

2446 

38.02 

l6o42 


1844 

52.21 

33,13 

2769 

44.0 

19.08 


2770 

51.83 

35.03 

3570 

46.42 

26.75 

4H 

3313 

73.23 

34.14 

4843 

43.2 

39.09 


4540_ 

79.55 

34.19 

5374 

43.9 

45.36 


62 06 

93.04 

33.68 

6429 

45.04 

59.36 


7619 

107.66 

34.84 

7089 

44.2 

72.81 


1626 

53.26 

32.0 

3 088 

34.8 

21.266 


2174 

58.16 

30.69 

3886 

35.9 

27.4 


2537 

62.06 

30.1 

4414 

37.2 

31.96 

2H 

3130 

66.19 

29.08 

5007 

37.9 

37.11 

3932 

73.35 

29.32 

5640 

39.2 

44.03 


43 42 

78.94 

30.98 

5836 

39.8 

47.96 


5604 

89.08 

30.71 

6630 

41.75 

58.37 


6372.5 

96.70 

31.19 

7038 

41.84 

65.51 


Table 4.14 : Multichip Natural Convection 

Open Closure 







— 

i 








u 

I O Q 








XI. 

n n 1 








■x-li 

Ki 





0 

B 


00 


ro m 

VO 

VO 

ro 

1 

m 

VO 

VO 

O CM 

VO 

in 

CM 


ro 

B 

If 

• 

o 

# 

00 


# i 

h- Ch 

® 

o 

# 

CO 

9 

LO 


If 





r4 

r-i 

r-i 


oo 

r- 

CO 

CM 



CO 

CM 

CM 

1 

xp 

r- 

VO 

CM Ch 

r- 

Ch 



CM 

0 

0 

• 

# « 

» 

♦ 

0 

<3 


a\ 

as 

CM 

l> CM 

CM 

CM 

r- 


I! 


r~f 

00 

^ VO 

(D 

O 

r-i 


Q 






r-i 

rH 


0 

B 

VO 

VO 


CM 00 

Ch 

Ch 


tH 

1 

00 

tH 

00 

00 O 

ri 

VO 

00 


tH 

« 

« 

» 

0 



9 

<1 

B 


l> 


Os VO 

O 

Ch 



II 



r-i 

rH CM 

00 

ro 




ro 

CO 

VO 







00 

!> 

ro 

3 

0 

o 

CO 


CM 


« 

» 

• 

9 « 

« 

9 

« 

03 


VO 

r-l 

rH 

r- r- 

00 

VO 

CM 



o 

CD 

Ch 

in CO 

00 

VO 

ro 



r-l 

CTv 

Ch 

VO r- 

Ch 

CM 

in 



CO 

in 

00 

rH ro in 

!> 

r- 






tM m 

r-i 

r-i 


CM 


00 

r“i 

CD 

CD CM 

o 


CO 

'::s 


LO 

ro 

Ch 

VO CM 

in 


m 



0 

« 

« 

« « 

• 

« 

9 



CM 



tn vD 

VO 

VO 



rH 

r-i 

r-i 

r-i T- 

r-i 

r^ 

r-l 

8 


CM 

00 

o 

2 

6 

Ch 

O 

VO 



« 

« 


# C 


« 

• 



o 

o 

o 

O o Ch 

O 

r-i 



CM 

CM 

CM 

CM CM rM 

CM 

CM 




ro 


ro ro 

VO 

VO 

ro 


rHI 

Ch 

VO 

CM CO 

m 

in 

00 



0 

0 

• 

» « 

6 

« 

9 



r-i 

00 


Ch 

o 

CO 

VO 



CM 

CM 

CM 

CM CM 

CO 

CO 

ro 



!> 

CO 

CM 



CO 

CM 

CM 


VO 

O 

VO 

4 

5 

6 

Ch 

O 

B 


« 

0 

• 

« « 

c 

@ 

9 



m 

a 

CM 

00 CM 

CM 

Ch 



CM 


10 

VO CO O 

CM 

ro 







r-i 

rH 

r-i 



VO 

VO 


CM ro Ch 

Ch 




m 


ro 

LiO VO 

O 

VO 

Ch 

B 


» 


0 

• 1 

« 

• 

9 

1 


r- 

rH 

G\ VO 

o 

Ch 

00 


i 

CM 

CM 

CO 

00 

tn 

tn 

VO 



r^ 









CM 

O 

00 







# 

o 

tH 

tn H 

m 

00 


CM 

1 

o 

0 

o 

0 

CO 

« « 

00 

e 

m 

9 

r-i 

o 


i 

CM 

o 

r- 

Ch 

CM 


CM 




tn 

00 

CM 

VO 

CO 

Ch 


i 




v-i tH 

CM 

CO 

CO 


0 

CNI 

II 

% 

X 

(n 

If 

P 


150 


I able 4.15 : Multichip Natural Convection (Closed Enclosure) 
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Table 4,18 : Multichip Hatural Convection 

Rise in temperature with heat flux 


AT 2 ZiT3 


48.9 

43.6 

1.3 

124.0 

122.1 

2.6 

227.4 

235.0 

4.97 

350.0 

359.1 

9.77 

594.0 

589.0 

15.22 

676.8 

685.9 

18.2 

1025.0 

980,0 

26.45 

1467.3 

1437.4 

31.85 

1762.4 

1768.2 

38.67 

2235.5 

2237.9 

49.75 

2992.3 

3008.3 

58.42 

3085.8 

3068.9 

62.72 


D = 4H, AR = 2 


2.72 

2.06 



6,55 

5,2 



12.32 

10.0 

1 


19,6 

16.2 

2 


30.15 

23.2 

3 

,D, 

34.95 

27,1 



51,0 

38.06 



64.9 

51.55 



78,2 

62.3 



97.27 

77.8 



116,92 

93,8 



122.27 

97.8 






Table 4® 19 1 Multichip Natural Coniroction 
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Table 4,20 


Multichip Natural Convection. 

Rise in temperature with heat fltax ^ 

3 


3H 



^2 

% 



^3 

42.8 

50.1 

53.81 

2.4 

2.53 

2.0 

187.8 

194.9 

194,85 

10.95 

12.72 

8.36 

459.1 

468.72 

464.6 

27,25 

27.4 

18,96 

752o7 

759.4 

729,6 

40,05 

40.17 

25.93 

1206,2 

1287.1 

1219,1 

60,6 

61.55 

42.73 

1272,1 

12 98*6 

1232,2 

64.22 

64.0 

42 ,36 

1855.3 

1873.2 

1808.4 

85.27 

86.2 

60.03 


D = 3H, AR = 2 
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Table 4o21; Single chixj Forced Convection 



39.36 

20.2 

18,0 

0 

3649,5 

17.45 

oo 

40,29 

19.9 

18.0 

250.7 

3161,0 

20.69 

0.050 

40,29 

19.6 

18.0 

368.1 

2669.7 

24.58 

0.019 

43,50 

19.25 

18.0 

452.2 

2092.8 

31.8 

0.010 

40,7 

19.15 

18.0 

553,3 

1927.2 

34.4 

6,3X10“^ 

40.5 

19.0 

18,0 

637.3 

1678.3 

39.5 

4.1x10“^ 

40,5 

18.83 

18.0 

721.5 

1395.2 

47.92 

2 •68x10"'^ 

40.7 

18.77 

13.0 

805.8 

1295.19 

51.41 

l*99xlo"*^ 

40.29 

18.76 

18.0 

889 .9 

1278.49 

52.08 

l^SlxlO**^ 

40.29 

18.70 

18,0 

973.9 

1178.24 

56.55 

1«24x10“^ 


= 40,64 w/m 


AR = 3, Average q 


Table 4,22 


D 


t Single Chip i^orced Convection 

II 


Vi 

■"c 

00 

Re 

Ra 

m 

Ra/Re^ 

67.12 

^5,2 

21.5 

0 

5724.3 

17.46 

cx> 

71 #64 

24.65 

22.0 

260,7 

4108.69 

26,03 

0,060 

6S .55 

25,15 

22,9 

438,1 

3452.7 

29.28 

0.017 

-3 

o4.31 

25,04 

23.4 

517.9 

2511.5 

37.06 

9.39r.lO 

“3 

63,27 

24.8 

23,0 

698.3 

2769.7 

36.99 

5,68x10 

-3 

71.32 

25-05 

23.3 

778.2 

2681.43 

39.42 

4,43x10 

-3 

6=.37 

25.2 

23. S 

85«.9 

2443.5 

41.02 

3.319x10 

-3 

71 .82 

2 6.35 

24.6 

939.4 

2627.56 

39.27 

2.97x10 

« ^ -.-3 

71.64 

2 5,6 

25.0 

1011.6 

2390.=? 

42.8 

2.33x10 


Average q = 68 .97 w/m 


AR 
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Table 4.2.< 


Single Chip Forced Convection 



104.2 

27.0 

21.5 

0 

8364.8 

18.1 

oo 

107.1 

25.35 

22.1 

243,13 

6479,7 

24.2 

0.10 

100.71 

26.2 

22.8 

436.9 

5168.9 

28.2 

0.027 

08*9 

25.15 

23.1 

614.47 

4630,5 

31.1 

0.0122 

101.8 

25.3 

23,5 

774.9 

4234,5 

34.8 

7 . 05 x 10 "^ 

99,62 

25.05 

23.6 

856,0 

3711.7 

33.9 

5 . 06 x 10 “^ 




Ta^le 4,24 ’• Single Chip Forced Cowectfon . 


q 

T 

c 

T 

oo 

Re 

Ra 

Hu 

R 3 ./R 0 

127 .0 

2 5 ® 5 

18,3 

0.0 

11329,4 

17,06 

CD 

125.15 

24,5 

19.3 

245.5 

8209.85 

23.27 

0.135 

125.12 

24.1 

19.3 

443.8 

7607.3 

25.42 

0.038 

125.12 

24,9 

20.1 

543.0 

7512ol 

25,36 

0.025 

129.78 

24,15 

20.1 

622.0 

6383.9 

30.97 

0,016 

329,0 

24.1 

20,4 

705.8 

5824.1 

33.68 

0.0117 

129.0 

24.63 

20.8 

785.1 

59 83.5 

32.49 

9.70x10 

^ ^-3 

129.36 

24.2 

20.8 

868.0 

5333.6 

35.88 

7.07x10 

-3 

124.42 

24,6 

21,1 

947.9 

5459.4 

34.28 

6.08x10 

^-3 

126.7 

24,4 

21.1 

1030.2 

5157 .27 

37.03 

4.86x10 


/ 2 

AR = 3, Average q = 126,96 w/m 
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Table d»2D : Single Chip Forced Convection 


Cl 

c 

oo 

Rg 

Ra 

Nu 

Ra/Re^ 

520.0 

44.95 

13.6 

0 

34517.4 

18.541 

OO 

534.0 

42.0 

19.6 

248 

29965.5 

22.46 

0.487 

520.0 

39.9 

19.8 

446 

27379.8 

24.44 

0.137 

529.0 

38.9 

20.3 

545 

25494.4 

26.88 

0.085 

529.0 

38.25 

20.2 

628 

24751.6 

27.89 

0.0627 

539.0 

38.3 

20.6 

710 

24351.5 

28.80 

0.048 

534.0 

37.95 

21.0 

792 

23338.9 

29.79 

0.037 

534.0 

37.7 

21.2 

875 

22744.4 

30.61 

0.0297 

549.0 

37.65 

21.6 

957 

22081.1 

32.33 

0.0241 

534.0 

37.3 

21.8 

1040 

21367.8 

32.57 

0.0197 


= 532.2 w/m^ 


AR - 3, Average q 


able 4,2 5 


Single Chip Forced Convection 


"”1 

T 

c 

oo 

Re 

Ra 

Nu 

Ra/ke^ 

3021.6 

53 .5 

19.0 

0 

49381.3 

20,92 

03 

1014.7 

58.25 

19.1 

248 

45377.4 

23.87 

0.7377 

1053.4 

55 #6 

19.5 

364 

43532.8 

26.45 

0.328 

1001. 0 

53.5 

20.0 

546 

40290,2 

27.67 

0.1351 

1001 .0 

53.4 

20.7 

628 

39202,0 

28.33 

0.099 

987.4 

51.65 

21.0 

710 

37258.0 

29.87 

0.0739 

994.0 

53,6 

23 .3 

792 

35713.2 

30.24 

0.0569 

947.0 

50.45 

22.8 

891 

33617.3 

31.73 

0,0423 

1031..5 

52.25 

24.1 

981 

33423.9 

33.93 

0.0347 

1001.0 

51.9 

24.6 

1040 

32421.2 

33.82 

0.0299 


/ 2 

AR = 3, Average q = 1006,66 w/m 



Ta’ole ‘i.n ; Single Chip Forced Convection 


q 


T 

oo 

Re 

1535.4 

36.5 

22,3 

0 

1511.2 

77.5 

22.7 

209.7 

1553.4 

75.25 

23.1 

379.5 

1519.0 

72.25 

24.0 

452,9 

1553.4 

71.7 

24.3 

537.5 

1596.3 

71.1 

24.3 

609,3 

1486.1 

57.5 

23.4 

688.2 

1613 .6 

58.7 

24.0 

760.9 

1553.4 

66,5 

24.6 

831.6 

1579.1 

66.9 

24.7 

873,4 


Ra 

Hu 

Rq./R© 

57371.1 

21.11 

oo 

52686.6 

24.61 

1.206 

50988.8 

26. 55 

0.354 

48151,0 

28.26 

0.235 

47445.0 

2 9. 42 

0.164 

47085.2 

2 0,64 

0.126 

45995,5 

30,46 

0.097 

45985.2 

32.55 

0.079 

43779.2 

33.5 

0,0633 

43916,5 

33.79 

0.0576 


AR 


3 , Average q 


1550.25 




oo 


Re 


1 } 79,3 

1337.3 
1 G 37 .3 
1847.0 

1810.3 

1319.3 

1337.3 

1837.8 

1337.8 



T 


33.2 

22.8 

0 

33.5 

23.4 

203.27 

3 1 .35 

23.7 

369.6 

80.1 

24.1 

525.7 

78.1 

24.4 

597.64 

7 5.5 

24.7 

569.49 

75.0 

24.7 

741.04 

75.4 

25.2 

310.2 

74.5 

25.3 

882,0 


61362.0 

22.45 

CX5 

56875.2 

24,80 

1.38 

54884.6 

26.62 

0.40 

52264.8 

29,28 

0.189 

50882.9 

30.14 

0.142 

49666.3 

31.30 

0.110 

49176,2 

32.22 

0.089 

48446,4 

32 , 66 

0.073 

47817.4 

33,36 

0 . 06 l 
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’able •:Sin 3 le Chip Forced Convection with a Protrusion 
on the Opposite Wall 


D 


q 

T 

Oo 

Re 

Ra 

Nu 

Ra/Re2 

103.43 

23,5 

22.5 

0 

8941 


OO 

108,5 

27.4 

23.1 

338.31 

6451 

24.16 

0,056 

109.9 

27,4 

23.1 

402.7 

6451 

24. 

0.039 

109.3 

27,35 

23.1 

499 ♦ 49 

6379 

24.76 

0.0256 

104,7 

25.6 

21.8 

569,10 

5845 

26,6 

0,0180 

103.02 

26,07 

22.5 

631,9 

5444 

27.71 

0.0136 

102.5 

25.3 

21.8 

683.5 

5399 

28.18 

0.0115 

100,8 

25.3 

21.9 

732.1 

5242 

28.62 

9 .78x10 

99,97 

25.27 

22.0 

780.7 

5040 

29.41 

8 . 28 x 10 


= 104.43 w/ra'^ 


AR = 3, Average q 


;b 


j-al Ic -1 

1 1 • 

Single Chip Forced Convection 

on the Opposite wall 

with a 

Protrusio) 

q 

*TI 

""C 

00 

Re 

Ra 

m 

Ra/Re^ 

127,9 

34.7 

22.5 

0,00 

17156.2 

17.71 

oo 

220^7 

35.25 

27.0 

325.8 

12480.3 

22.42 

0.118 

224.6 

35.6 

26.7 

389.0 

12102.0 

23.75 

0.079 

219.15 

35.9 

27.8 

480.3 

10911.7 

25.41 

0.047 

220,44 

35.9 

28,3 

541,6 

10207.7 

27 .22 

0.034 

219,8 

35.9 

28.4 

603.2 

10067.3 

27.50 

0,0276 

22 0,93 

35.3 

28.2 

651.2 

9594,7 

30.03 

0.0226 

219.8 

35.35 

29,0 

695.0 

9166.2 

30.09 

0.0189 

232.18 

36.5 

29.9 

738.0 

8732,2 

32.91 

0.0160 


/ 2 

= 223.5 W/m 


U 


AR = 3, Average q 
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*1 «* 

Single 

on the 

Chip Forced Convection 

Opposite Wall 

V7ith a 

Protrusion 

a 

m 

1.C 

oo 

Re 

Ra 

Nu 

2 

Ra/Re 

763. 1 

63 . 3 

28.2 

0.0 

36939.7 

19.77 

oo 

763.41 

59.1 

29.1 

303.4 

32578.53 

23.24 

0.354 

780.47 

57.9 

29.8 

361.87 

30713.55 

25.21 

0.235 

768.41 

55.5 

29.9 

451.6 

28559.42 

27.34 

0.140 

7 68 .41 

55.0 

30.4 

509.8 

27485.65 

28.45 

0.106 

774.43 

54.85 

30.7 

567.9 

26971.98 

29.20 

0.083 

780.47 

54.6 

30.9 

611.67 

26497.20 

29.99 

0.070 

861.10 

54.4 

31.0 

699.2 

26192.8 

33.52 

0.0536 
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f ’1 

4 

.32; 

liUltichip Forced 

Convection 





"3 

Re 

Rai 


Ras 

NU 3 


133.4 

174.2 

00 

2865 

11.25 

3083 

9.59 




134.5 

170.3 

142 

2414 

13.28 

2555 

11.52 

1 

i 

3 

H 

187.0 

157.3 

243 

1971 

16.63 

1855 

13,88 

2 

3 

4 

182.4 

lo3.9 

348 

1826 

17.61 

1942 

14.81 

3 

3 

H 

178.2 

15 9.5 

448 

1647 

19.17 

1794 

15.66 


44 


190.5 

1 oS . 4 

524 

1526 

22.17 

1466 

20.44 


11 



AR = 2 

2 

Average q- = 184,9 w/m 

2 

Average ” iSo.S w/m 


Table 4.3 3 ; Multichip Forced Convection 


^^l 

% 

Re 


Nu^ 

Ra3 

NU 3 




606.4 

586.3 

00 

5818 

14.83 

6358 

12.77 


— 1 


620.5 

588.8 

135 

5597 

16.16 

5993 

14.04 

i 

J 

H 

636.2 

592.4 

234 

5000 

19.04 

5396 

16.13 

1 

□ 

c 

652.02 

625.5 

337 

4925 

19.97 

5215 

17.85 

3 

3 

H 

633 ,6 

614.3 

432 

4795 

20.09 

4966 

18.66 


11 


620.6 

608.2 

507 

4710 

20.15 

4780 

19.40 




AR = 2 

2 

Average q- = 627.8 W/m 
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Table 4 : Multichip Forced Convection 


qi 

'^3 

Re 

Rai 

Nu, 

Ra3 

NU 3 

1323.3 

125:. 6 

0 

9430 

15.74 

9871 

13.68 

1333.7 

1361 ,4 

127 

8760 

17.87 

9093 

15.84 

1435.4 

1 3 60.3 

226 

9038 

19.58 

9545 

16.96 

1545.0 

1337.1 

313 

8902 

20.05 

8917 

18.13 

1525.2 

1415.8 

406 

8520 

21.23 

8859 

18,51 

1425.6 

1305.8 

480 

8349 

20.48 

8399 

18.59 


AR = 2 

2 

Average q, = 1432,8 w/m 

A 2 

Average q 2 1331,6 w/m 


1 □ 

2 □ 
3 


H 

C 

H 


Table 4.35 


: Multichip Forced Convection 


**1 

^3 

Re 

Rai 


Ra^ 

NU 3 

2380.0 

2380.1 

0 

12605 

16,40 

12977 

15.07 

23 J4.7 

2247.0 

120 

11416 

19.26 

11771 

17.50 

2392.7 

2417.0 

209 

11662 

19.05 

11835 

18.59 

2279.8 

2223.3 

300 

10973 

19.73 

11138 

18.63 

2418.1 

2331,3 

385 

10860 

21.05 

11001 

19.74 

2469.5 

2307.09 

453 

10765 

22.04 

10795 

20.74 

AR " = : 

2 






Average 

<^1 “ 

2373 

.0 w/m*^ 

2 




Average 

qS = 

2317. 

5 w/m 





Multichip forced convection® Temperature varxation '‘rith Reynolds number 



Table 4.37 * Multichip forced convection. Temperature variation \rith Reynolds number 




Multichip forced convection. Tempera-cure variation i^/ich Reynolds numbe: 
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Table 4.33 . Multichip forced convection. Temperature varratron uith 

Revnolds number 
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Table 4.41 * Multichip Forced Convection 
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Table 4.A2- * Multichip Forced Convection 
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Table 4.46 s Multichip Forced Convection 
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Table A-A8 s Multxchip forced convection. Temperature variation with Reynolds number 
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Tcble A‘5X : Multichip forced convection 
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CONCUJSIQK 


TTip prcstiHt work concerns the study of free, mixed and 
forced convective h'-at transfer from the single discrete heater 
and ar. errny three heaters in the vertical configuration. 
Rerultr, obtained experimentally show the follow; 

.^^’Inple . Me a ter; ?lov.’ patterns in the vicinity of a single 
heater are complex and a universal self similar relationship 
be twee:-, the average heat flux and the wall tenperature is not 
possible, '.’u-.rclt numbers for a protruded heater generally 
exceed values quoted for flush heater. As the heater length 
reduces the average h'usselt number increases owing to the 
approach tov.'ard.s point source. Heat transfer is unaffected 
by the presence of a neighbouring wall and infact goes down 
v.'hen the flov.* path is baffled. Unlike the problem of a flush 
heater buoyancy eff-^cts do not enhance heat transfer nKonotoni- 
cally ove-r and above the forced convection value. Adverse 
stratification effects and the stabilizing effect of buoyancy 
on recirculation can lower the Wusselt number over portions of 
mixed convection regime. 

(2) Array of Three Heaters : The heater placed at the lowest 
ixosition shows the best thermal performance while the next to 
it shows the worst performance. Subsequent heaters do not show 
further deterioration mainly because acceleration of flow and 



disturbances introduced by the heaters themselves. Consider- 
able advantage is gained when an unpovjered heater is placed 
between tv^o powered conponents. As in the problem of a single 
heater the role of the neighbouring passive wall is minimum. 

In natural convection the plume is diverted towards the 
passive wall and this suggests the use of cooling strategies 
in engineering practice. Experiments for an array of heaters 
shov7 that a definite heat transfer augnentation is possible in 
mixed convection, Hoxvever the thickening of the thermal boundary- 
layer effectively reduces this augmentation as one moves verti- 
cally upwards. 
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1 

k - '.C AL [ALYS IS 

Bias and random errors are the t\7o types of errors 
ansociatod 'ith experimentation* Bias errors are relateo ro 
tns raeasuring probes and measursLaent system. With the most 
raocern systems having a self calibrating capacity, these errors 
can loe minimized. The random errors are difficult to control 
sj.nce they arise from a number of factors that simultaneously 
afreet the experiment. These include room temperature changes 
supply voltage fluctuations, air currents, building vibrations 
and roughness of nominally smooth surfaces. These errors 
inherently e::ist in any experiment rendering the experimental 
data as fundamentally irreproducable . For a large number of 
experi.aen tal runs, provided the experiments are of good qualit 
results are distributed about a stable mean value. Such a 
distribution is called scatter and all efforts are needed to 
maintain a low scatter value. Overall improvement in the 
quality of measurement systems and test cell along ■with 
control of external factors influencing the experiments may 
result in a reduced scatter. The level of scatter is determi. 
by a series of experiments dedicated to determination of scat^ 
only. These experiments involve performing runs in identical 
conditions on different days, runs with measurement systems 
having identical specifications and runs on different apparat 
that look alike. 

Error Propagation 

Due to scatter an uncertainty is always associated 


With a measured variable. If a new quantity is defined on 
the basis of this variable then the uncertainly involved 


may also affect the value of the new quantity. If 6 
the uncertainty involved with a variable and y be the 

function of x. then uncertainty in y, 6 ^ win le 

X V 


C 2 

i 


ax. ^ 

X 


1/2 


On the basis of error propagation formula and a series oj- 
experiments dedicated to the determination of scatter, 
c recent the level of scatter for this v7ork. 


TJatural Convection Experiments 
we have , 

2 

hL _ ql» - V h. 

NU = ^ ~ 

V is the voltage applied, volts 

L is the characteristic length of heater, m 

2 

A is the area, m 

R is the resistance. Ohms 

k is the fluid thermal conductivity, W/rnK 

iir is the temperature difference between heater surface 

and ambience, C • 

For this problem the resistance R of nlchrome wire with 
length L and area A of the heater can be considered as 
constant. The uncertainty in the Husselt number will 
depend upon the uncertainty in the measurements of vol*^ 9 
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■cemperature and thermal conductivity of fluid. By error 
propagation formula we get the uncertainty 6 Hu in Husselt 
number as 


^ gq 
Hu 


[ 4 (-^)^ + (-4#)^ + (4^)^ ] 


T 1/2 


AT 


k 


and uncertainty in Rayleigh number as 


^ = r ^1/2 

Ra ^ 6 AT ^ rf -• 


V « 

The prefix 6 indicates the maximum variation of any 'niantity 
from its mean value. The denominator quantities are the 
mean values of different variables. We have carried out 
the exoeriments with nearly same input heat flxjx on different 
days in identical external conditions. Table All presents 
the data for these exoeriments. From table we calculate - 


V = 28,58 V 

mean 


11.84 Watts 


.o. 


mean 


T , . = 62.03 C 

m(mean) 

AT ^ = 61.92 

mean 


.-3 


. = 28.83x10 " w/ltlK 

Kraean 

Maximum variation in - 

voltage fron its mean value - 0.21 volts 
rff from its mean value: -2.32°C 


k from its mean value: -O.SXfO""^ uncertainty 


in the Nusselt number can be calculated as. 
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^Nu _ 
Nu 


E4(' 


0 . 21 >2 

28.58 




f-2 ,«.32 2 

'■ 61 . 92 ^ 


0.3 v2 ni/2 
28.83^ J 


£Nu 

llu 


= 0.041 


The uncertainty in the Nusselt number is + 4.1%. 

Uncertainty in Rayleigh Humber 

3 


Ra 


B A T g L' 


va 


Value of B depends on the heated surface temperature 


CT^ + 273) 




The value of g and L can be consicered as constant, 


-»5 2 

a = 2,77 X 10 m /sec 

mean 


4.T = 61.92 C 

mean 


-3 


mean 


y 


mean 


= 2,728 X 10 /K 

-6 2 . 

18.95x10 m /sec 


Maximum variation in 

3 

fl from its mean value s 0,022 x 1^^ A? 


-5 2 


a 


from its mean value : 0.1x10 m /sec 


5 ; from its mean value 


6 2 

0.16x10" m /sec 


isT from its mean value: -2.32 C 


How, 



6 Ra 
Ra 



0.022 x 2 

2,728^ 




(^2 2 
^ 61 . 92 ^ 


48.95 


+ 


'' 2.77 




0.053 


The uncertainty involved in Rayleigh nunber is 5.3%. 

2 

We have Q = V /R 
= 6 Q/Q = 

Maximum variation in voltage is 0.7?^ and the corresponding 
variation in heat input is 1.4%. 

It can be stated that for a variation in supply voltj 
within + 0*7%. 

Power input varies within + 1.4% 

Nusselt number varies within + 4.1% 
and Rayleigh number varies within + 5.3%. 

We have also carried out experiments to determine scatter at 
different power input the results of which are stated below. 

Power Input Q = 5.6il Watts 

For a variation in supply voltage within + 2.09?^ 

Variation in power input Q is within + 3.56% 

Variation in the Nusselt number is within + 5.6^o 
Variation in the Rayleigh number is vjithin +8.3% 

Power Input Q = 2.801 Watts 

For a variation in supply voltage within ± 1.735% 

Variation in power Q is within ± 3.49% 

Variation in the Nusselt number is \’ithin ± 7.25% 

Variation in the Rayleigh number is within +6.87% 



We observe that the uncertainty in the Nusselt and the Rayleigh 
number increases with a lox-^ering in input po'v7er. Tn.i=, is 
because of the fact that as the power input is lo\ 7 ered the 
natural convective flow gets weaker and weaker and are likely 
to be affected more by the external influences. For the present 
w^orh the overall uncertainty in the Nusselt number can be taken 
as 5% and that in the Rayleigh number as 6%. 


Forced Convection Experiments 

A similar analysis has been done for the forced conve- 
ction e> 5 )eriments. The uncertainty in the Nusselt number and 
the Rayleigh number is found to be 4%, The uncertainty in 
Reynolds number varies \-d.th velocity. At lower velocities 
(U^ 0.3 m/s) the uncertainty involved is high and is about 
6%, At higher velocities (U^O.65 m/s) the uncertainty 
involved is low and is about 3%. 
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WORKING PRINCIPIiS OF INTERFEROMETER 

The interferometer works on the principle of inter- 
ference of light. Interference is the phenomena of redis- 
tribution of light energy which takes place when two light 
waves originating from the same monochromatic source get 
superimposed. Ijight is characterised as the visible spe- 
ctrum of electromagnetic waves. For a monochromatic source 
only one wavelength is significant and the propagation of 
light originating from a monochromatic source can be expre- 
ssed as 


E 


A Sin 2%{vt - 



where E is the electric field. A is amplitude v is 
frequency# x is wavelength# t is time and x is distance. 

The superimposition of two v/aves having a phase 
difference 0 from the same monochromatic source can be 
expressed as 

+ E2 = A [ Sin 2 Ti 3/ t - -^) + sin( 23i(tj;- -f~) + <SS ] 


The amplitude of both the waves is same because 

they are produced by the same source. The sum of the 


electric field vectors is# 


Ei+ E2 


Sin( 2 Tivt 



2 A Cos 



The intensity of a light wave is proportional to square of 

2 2 

amplitude. The intensity of combined beam is 4A Cos ^ . 
This can be plotted as 



It is clear from the plot that at some points there is an 
addition of intensities which will form the bright spots 
and at some points the intensity is zero which will form 
the dark spots. This redistribution of energy leads to 
an interference pattern in the form of dark and bright 
regions A and B respectively . The spacing between two 
consecutive points of maximum or minimum intensity is 
termed as fringe spacing or fringe width. 

An interference pattern is stable only if the 
conditions of spatial and temporal coherence are satisfied. 
For spatial coherence, each wave should have traversed the 
same length x and for temporal coherence the phase at a 
point and the phase difference between two points should be 
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constant with time. 

No light source is perfectly coherent in all respects. 
The degree of coherence is judged by its coherence length 
which is the distance over which the waves emitted by the 
light source satisfy the conditions of coherence. Laser is 
a source which gives the coherence length of the order of 
few kilometers. 


In the present study we have used a He-Ne laser as 
a light source . It has a continuous wave output of 35 mVI 
and a wavelength of 632.8 nm. Figure A 2 2. shows sche- 
matic diagram of interferometer. 



E>£j3Tyi SplitttT 


Schematic diagre oZ 'C'.-Ihrn :r ro <'--r 

In the Mach-Zhender interferometer, the interference 
pattern is formed due to a phase difference between the test 
and reference beam. 

The index of refraction of a homogeneous medium is 
a function of the thermodynamic state, often only density. 
According to Lorenz-Lorenz relation, the Index of refraction 



of a homogeneous medium can be obtained from 


1 

P 


2 

n 

n 


- 1 

+ 2 


When n ~ 1, this reduces to 


constant 


n-1 


0 


c 


=> 


= constant 

d P 


For small changes in temperature (upto 20°C in 
air) density and temperature are linearly related. 


= constant can then be \:ritten as ^ = constant 
and any change in T will have a simultaneous effect on n. 

If a beam of light passes through a medium of varying temper- 
ature then the optical path length traversed by beam will be 
more than the geometrical path length. 



L 

Path length PL = / n dz > L since n> 1. 

o 

Let the test beam moves through a region of variable 
n and reference beam moves through a region of constant n 


then 
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L 


PL = 

PLi 

- PL 2 = 

/ 

0 

(ni “ n2)dz 

= 

dn 

df 

L 

/ 

0 

( ei - 

f2> 

dz 

= 

dn 

dT 

y 

0 

(Ti - 


dz 

If T^ 

1 

is a 

two 

dimensional 

, temperature field then 


tne integral simplifies to 


PL = (Tj - T2) 


dn 

dT 


. L 


A path length difference of A generates one fringe 
shift, the At for one fringe shift can be evaluated as 


At 


For air (f") 

uT 


\/ (I 


dn 

dT 


) 


0.927 X 10“V°C . 


For He-Ne laser \= 632.8 nm. Hence LAT par fringe 

shift for a He-Ne laser is 0.632 m°C in air. The temper- 
ature change per fringe shift depends upon the length L of 
test section. More the length, more sensitive the inter- 
ferometer will be. But large values of L lead to high 
refraction errors. Each fringe is a line of constant phase, 
refractive index, density and temperature and is an isotherm. 

Consider a ray travelling in a medium of varying 
refractive index and temperature, 

t I ^ 

^ “H- Fig.A.2.4 



The effect of refraction is to displace the beam towards 
regions of higher refractive index n# hence higher density 
and loviBr temperatures. This increases the optical path 
traversed by the beam. 

Interference patterns now formed are due to the 
combined effect of difference between the average refractive 
indexes for the test and reference beam and an additional 
path difference introduced by gradients n in test cell. 

The fringe shift e can be considered to be. 


%otal '^true '' refraction 


is the fringe shift caused due to the differ- 

^true 

ence in average refractive index. 

£^ef is fringe shift introduced by the refra- 

ction of test beam when it passes through test section. 


Temperature is related to the difference in average 
refractive index for the two beams alone hence any fringe 
shift introduced by the refraction is a source of error. 

mup, ^ can be calculated by the formula 

^refraction 


^ref 



( 011)2 

MT 




2 


AS the length L increases this term increases as I, . In 
practical values of I, > Im the refraction errors are subs 


tantial 



lya 


In the present study the maximum length of test 
section is 0.2m. Hence the magnitude of refraction error 
is small. To find out the local wall heat flux the fringe 
spacings closest to the wall is measured and heat flux is 
given by 


q(x) = k(fluid) 


AT/fring shift 

e 


and 


Nu 


AT/frinqe shift 
e Cx) 



oo 


) 


When the heat flxix is large the wedge fringe mode is preferred 
as the initial setting of interferometer. In this mode the 
fringes are initially made perpendicular to the surface to be 
heated. 



Pig.A.2.5: Initial 




As the temperature of the surface rises the fringes 
get shifted close to wall. 



Fig.A.2.6 : Wedge fringes at heated state 


This fringe shift is due to an increase in the tem- 
perature and is a measure of temperature change in that 
region. The curve A'B‘ represents the temperature profile 
at an x location A* and the slope of this curve at wall is 
proportional to the wall heat flux. 

This slope can be calculated from a photographic film 
and the heat flux calculated as 

a = -i. 

^ K 3yiy=0 

A curve can also bo fitted through some points away from 
from the heated surface and approximate temperature profile 
be obtained. This method has been elaborated in Chapter 3. 


can 
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Appendix 2 

CALCULATION OF RADIATIVE HEAT LOSS FROM A SINGLE CHIP 

HEATER 


Here we estimate the radiation heat loss in a black 


enclosure. The size and geometry of the enclosure are given 



ux 


Heat exchange between the aluminium strip and walls 
of the enclosure can be found out using the formula. 


Q 


c-e 




1 -e. 


1 -i 


’c- e 


v/here subscript: c refers to the heated strip and e refers to 
enclosure . 

£, is emissivity and A is area the shape 

factor between the chip and the enclosure. This formula 
assumes that the temperature of surfaces is uniform. In 
the present case the temperature inside the enclosure 
varies with position. To circumvent this problem the encl- 
osure is divided into many parts each of which is assumed to 
be isothermal. The heat exchange between the strip and each 
isothermal part is determined separately. The shape factors 
for all the geometries are determined using standard tables. 
The enclosure is divided as shownbelow. The associated shape 
factors are given along with the figures. 



^ 2-12 ~ ^ 2-14 

^ 2-13 


0.196 


( 1 ) 


0.5 
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( 6 ) 


fo = 0.0494 

2-ambxent 

Fig, A. 3. 2 

Heat exchange in enclosure = Hear exchange in the configu- 
rations (1+ 2+ 3+ 3+5+6) , 

for the polished aluminium e* strip can be taken as 0*2 
(from the standard books) • 

for the black enclosure is assumed as 1 because it gives 
•che maximum radiation loss* 

In the present work. 


^ 3 


= 0.2 


Table gives the temperatures at various positions 

in the enclosure and the amount of heat lost from the strip 
purely by radiation. 

Sample Calculation 

The heat lost from surface 2 to surface 12 can be 
calculated as shown belowJ 


Q 


^ 2 ^ ^ ^2 


2-12 


1 -e, 


1 -e 


12 


" 2-12 


12 


^2 

^12 
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^2 - 2«5 X 18 X lO”*^ 

^ 2 - 12 = 0 - 19 ^ 

= Planck's constant = 5,66xl0“^ — 

= 95,2°C = 368. 2Aj T^2 52°C = 325A 

2.5 X 18 X lO"^ X 5.66 x lO"® (368.2^-325"^) 

Q2-12 = — 

1-0.2 _J 

0.2 ^ 0.195 

^2-l2 ~ 0.202 watts. 

Similarly heat exchanges for other pairs of surfaces can 
be calculated. 

Experiments have also been performed in a \;hite 
enclosure to determine the role played by radiation in 
chip cooling. This is discussed in Chapter 4. Results 
presented in this work have been obtained in a white 
enclosure. 
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Radiation heat loss for a three dimensional heater. 
Emissivity of copper block can be assumed to 


e= 


Shape factor 


g' = 5.66 X 10 w/m^ k"^ 


= 53.2 C = 53.2 + 273 = 326.2 A 

w 


^oo = ^ ^ 


%ad 




= 0.15 X 5.67x10"® (326.2^ - 305^) 


Iraa = 22.78 w/n,' 


^in = 1S26 w/m"^ 


Irad = 1-^°° %n- 


It is seen that amount of radiation loss is very small in case 
of three dimensional heater when when compared with a two dimen- 
sional heater. This is due to the smaller surface area of three- 
dimensional heater and lower maximum temperature reached for same 
heat input because of higher heat transfer coefficients* 
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Append i.;c 4 

l^TUMJRIGAL CALCULATIOH OF HEAP TRMISS'SR FROM A 7 LUSH 

HEATER IITRIDB A GHAHI^TEL 


Here xve present mathematical modelling of a flush heater 
mounted on one v/all of a channel. The other v;all is assumed to 
be insulated and the top portion is modelled as an isothermal 
reservoir. The configuration is shoxm below* 

X.- l 

1 

X xchl. 

y x=o ^^‘■7 

Eia. A. 4.1 : Gonf iguration of flow and heat transfer 

ov6r a flush heater 

The channel is assumed to have infinite dimension 
in the third direction. Therefore the flov? is t\jo dimen- 
sional. A flat velocity profile has been taken but the 
scheme can easily be modified for a parabolic velocity 

profile • 
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orTiul-tion 

j-’he straaj! function 'j# vorticity u and temperature T 


the 

governing differential equation for a 

2 -D floxv. 


Toe equations are given as 



t + '? = - u 

XX yy 

(A. 4 ^ 1 ) 

a. j. 

J- 2 

-i-Uei' + v«fA~ 6 j +W + !;r“ H.e 

U.a.j^T V.U!^ yy Pj. 

( A # 4 # 2 ) 


t u T^ + V Ty = (T,^ t Tyy) 

(a « 4 ® 3 ) 


\ i tZ‘ ^ 


u = and V = - f • 

y 


jounaary coow'itions 

y SS 0 0 < X < hi/ 

n > hi, 

0 and Y “ 1 

0 and T = 0 


T = 1 


T = 0 

y 


y » edgey , 


0, '? = y « 


0 , (0 =t 0 and 


:z = edgex/fj^ 

In the present study edrjex = 5 , edgey - 1 and hi 1 * 

The method of false transient Is used to calculate 
steady state flov, and temperature patterns and average 
^.usselt number. In the n.th^ of false transient eguatrons 
(^.4.1) to u.4.3) are cast in the following form 


J M 


3t 




J 


(A. 4. 4) 

por incompresslbln fl» u^ f v^ = 0 we have 
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+ (vj2f) 

V 

u + ■v0^ 


In G ruction (A, 4, 4) 

0 

stands for 

any of the dependent 

\ C. 1 jc>i <.'0 and a,b,c cepend 

on 

the choice 

of 0. This corres- 

'on ence is given m Table 

jra 

below: 


^ a 

b 

c 

s 

1 

0 

1 

(l) 

1 

1 

1 

Ra . 

' “T 




Pr Re 

T l/Pr 

1 

1 

0 


Above 


nx UG 


•fforcntial equation is solved using numerical tech- 


bj.f£crence Approximation 

finite difference approximation for differeTcial 
•'■oa-^cjLon is given below. 

Present scheme uses 

a. Pull differencing for the transient term for 
^h+1 ^ 




at 


b. Central differencing for diffusion term 




Piti - ^1-1 - 


AX 


ay 


c. Upwinding for the advective term: 


V. u i2f = 


"r ~ ^1 ’"l 

Sc 


- Vb 


ay 
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i.‘- uh- carrent time level at which bhe solution must 
. .r. j-nc .» "ite assume constant step sizes , iit, ax and 
u.^ion mA) is no’j written in the finite diffcr- 


'c _ - u zr. follo\’SJ 


/ 




n-- : 


/ 


Zr 

a r V Zj 


nf 1 


h^[v . ujf] 


(Jf"-!) 


{ A ® ’4' # !d } 


in 




he? 




- vcctivG terms, r,l,u and v refer to right, left:, 
oclo-- reference point (i,j) respectively with respect 
local floiT direccion. This is shov/n in Fig. A. 4, 2 



jc onrcriiX a second order upx.dnding method for the advective 
btr .0 iv , u, /) below. The upvdnding involves backward 

(f»a 

‘ilferencing with respect to local flow direction. Define, 


u 

r 

11 

(u « 

“i+ih ' 



_ 1 
J 

(Ui 




„ 1 

- J 

(Vj I- 



p 

_ 1 

2 

(Vj + 




= 

efj, if 

Uj. a 0 , 

- d if u < 0 

“ ^i +1 ^ 



^ 1-1 

if u 2 ^ 0 ^ = 

Zi if ^ 4 . ^ 




V 0 # ^ 

of, ^ if V < 0 and 
^j +1 u 



:nni 


c + 


2 lit 

( A x) ^ 


2 /i i: a' 

(/Li y) ' 


0 


at b 


— -r - -2 V”' + ^ 


J2( 




Uc.fcxon (>k«4,o) is cast as a matrix equation the identical 
-orr. o. this equation for each dependent variable (0 = w, ^ and 
-/ o-or some simplification in computer programming 

oriort. In the method of false transient c^ = 1 for the 

0 

^-rvau 'unction and C = 0 when the real transient solution 

lUst oe obtained. Except when changes in the cIott quantities 

ic ra ''id the approach to steady state can be accelerated by 

0 

cr.oosing c >>1, The precise value of this quantity must 
j^e decided based on numerical experimentation. The choice of 
c > 1 IS equivalent to the successive over-relaxation (SOR) 
inethoa that is used with the Gauss-Seidel scheme for matrix 
inversion. 

The equation (A. 4, 6) is repeatedly applied for each 
node and for eecli time step. Further, within a time step 
iterations are required for matrix inversion (in the Gauss- 
Seidel sense) as well as for the non-linearity of the govern- 
ing equations (Equation A.4,2) . T'/pically 5-10 iterations 
are required per time step when the solution at the nth time 
level is used as the initial guess for the (n+l)th time level. 

A convergence criterion of 0.1 — 0.01 percent change per 
iteration in any of the flow variables at any node is generally 


used. 
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The grid used for a particular problem is chosen on 
the basis of several factors. Since a highly refined grid 
can lead to a high cost of computation the grid must be coarse 
but, a, must be fine enough to resolve boundary-layers, and 
b. must give a small excess energy balance at steady state. 

In dimensionless form, heat transfer rate is given by 

the Nusselt number, Nu. Nu is conventionally defined as hL/ 

k(fluid), where L is a lengrh scale and h is the heat transfer 

coefficient. For this problem 13u reduces to (- |^jy_o^ 

the average value of I^u is, 

fiL 

Nu = / Nu dx. The Nusselt number is specified as 

o 

- *ju (Rs, Pr) . 


Grid Refinement 

In wall-driven flows it is necessary to discretize 
the wall region with a sufficiently large number of nodes to 
resolve the boundary- 1 ayers . If a uniform grid is used a 
large number of nodes would be present in regions where velo- 
city as well as temperature vary gradually in space. This 
computationally uneconomical. Hence it is desirable to refine 
the grid locally near the wall while the central core has 
only a few nodes in it. This is schematically shown in Figure 







(k / 
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rciinenent referred here keeps the total nuober of nodes 
’.oo. ang>'d ■'mile it is selectively redistributed over the flov? 
-i0..tain. On a grid where the spacing ay is a variable conven- 
tional lor'^iulao for derivatives given earlier do not hold. 

M grid transf orraation approach is described to circumvent 
cnic problem. 


lor the channel considered here the boundaries y - 0 

... to :i = and y =«4jt>for all x are respectively heated and 

the y-direction in 

cooled. Hence we consider grid relin^- inA?hich the grid 
opoCing IS non-uniform. Let t) be a transformed plane in which 
the arid aopears uniform. Conventional finite difference for- 
-nuiac are applicable in this domain. Without proof we state 
tnat the following equation accomplishes grid refinement in 
the y-plane for uniform grid in the h -plane. 


dY 

dr\ 


Q (■ 




dri 


= 0 


(A. 4.7a) 


..-here qM = - A sign ( tl- \) ex,(-c | n-Hi ) 

i=l 

1 = 1 is the lower wall, 1 = 2 Is the top wall, hi = 0. 

h = 1 and A and C are constants. Typical values of these 

^ - n 0 At 25, Equation (A. 4. 7a) must 

constants are, c = 0.2 and o < A< . u 

be solved subject to the boundary conlitions, 

y ( 0) =0, y (h = 1) = 1 U.i.lh) 

This solution can be accomplished using finite differ- 
ences and is not described here. 




J 
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Equation (A.4«7a) is equivalent to, 

2 

d T) _ 

- Q when y is used as an independent variable. In 
dy 

two dimensions grid is generated using equations of the form 
2 _ 

yy called as elliptic grid gener- 

ation and is a seoarate field of study. Use of differential 
equations of the form given above guarantees that the movement 
of the grid lines within the physical plane is gradual as \\rell 
as confined vTxthin the true physical domain. 

The effect of grid refinement is shown in Table A. 4, 2 

below. 

Table A. 4. 2 Grid Refinement. C = 0.2 

ATI = 0.1 All = 0.01 


A = 0 

Ay (min) = 0.1 

mniiiQQQimQi 

10 

0,082 


25 

0,065 

0.00617 


Viewed through the coordinate T| the governing equ 
at ion (a. 4. 2) gets transformed to. 


Q (Tl ) + >) jZf r]ril *’ 

(a J + (v j2f ) / y 1 + (A. 4. 8) 

L I T) T]-* 

Here I = x and we use the transformation rules Yt^ 

and jZfyy = where ^ = ( l/y^ )^* wall tem- 
perature gradient is evaluated as 9T/®yjy_Q = T^/ y-rilT]:; q* 

The discretized version of Equation (A. 4. 8) is similar to 
Equation (A. 4. 6) and given below. 
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At a' 


0 


’“ 2Ay 




0 


m-1 


j+1 "j-ri 


At b‘ 


At b' 


0 


0 


AX 


/i t a' 

(a y)' 


'^3 “r - >' 


“>3 ''u - *=4 ’'b>“'^ J/ °®“"’ 


rri-1 


(v ) . A.y 

^Tl^J 


(A*4 .9) 


Denm 


0 

c + 


J 




2 A t a 2 A t a 

(Ax) (Ay) 

.0 


a t b 


aX 




)?( 

1 at b /, •u. .. \ n4" 1 


(y ) . ^y 

T] j 


(bi - bj v^) 


In Eq, (a. 4. 9), ax is the grid spacing in I ( = x) direction# 

Ay is the uniform grid spacing in t) direction and j is an 
index in the direction. 

It is noted that the grid generated from Equations 
(a. 4. 7) using a finite difference method will itself be affe- 
cted by discretization and round-off errors. The flow solution 
obtained from such a grid can be in considerable error since 
the metric derivatives such as dy/d'^ are involved in the 
flow calculation. Hence one must solve for y(T)) as accurately 
as possible. Usually y is determined on a grid with twice as 
many points as what is required for flow calculations. 
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Three programs are to be run separately in a fixed 
order to arrive at the solution 

1 For generating grid. 

2 For initializing the flow and s temperature field. 

3 For solving the flow and temperature equations using 
the method of false transient (c = 1) . 

The results are obtained by varying the values of Ra 
and -i^e in rhe program each time* The solution at a lower 
value of ^^a can act as the initial condition for the new 
problem* For the present case the runs have been made 
V 7 ith^ a grid size of 41x11, with c = 0.2 and A = 50. Conver- 
gence criterion of 0*1^ is used, Lt is 0.01 and time step 
marched for each case is 800, 



ii3r'4 





